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y  This  paper  describes  some  experimental  observations  of  free-surface  flows 
arising  when  fluid  poured  over  the  end  of  a  flat  plate  into  a  reservoir  well 
below  the  plate.  This  class  of  flows  was  found  to  be  highly  unstable  with 
over  seven  hundred  qualitatively  different  flows  being  observed  in  the  experi¬ 
ment.  At  very  small  values  of  the  flux  the  liquid  fell  from  the  plate  into 
the  reservoir  in  the  form  of  droplets  which  developed  periodically  at  a  number 
of  sites  on  the  underside  of  the  plate.  As  the  flux  was  increased  these 
'sites'  were  able  to  sustain  continuous,  unbroken  streams  and  combinations  of 
drip  sites  and  continuous  streams  were  possible.  At  still  larger  values  of 
the  flux  there  was  a  tendency  for  the  liquid  to  fall  in  'sheets'  and  several 
combinations  of  'sheet  fi  ws'  and  continuous  streams  were  observed,  some  of 
which  flows  exhibited  a  rather  chaotic  temporal  behaviour.  At  even  larger 
flux  values  the  flow  resembled  that  of  the  classic  waterfall,  but  here  there 
were  also  some  unexpected  instabilities  where  the  attachment  line  of  the  free 
surface  of  the  liquid  with  the  plate  developed  corrugations. 
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ABSTRACT  (continued) 


and  continuous  streams  were  possible.  At  still  larger  values  of  the  flux 
there  was  a  tendency  for  the  liquid  to  fall  in  'sheets'  and  several  combina¬ 
tions  of  'sheet  flows'  and  continuous  streams  were  observed,  some  of  which 
flows  exhibited  a  rather  chaotic  temporal  behaviour.  At  even  larger  flux 
values  the  flow  resembled  that  of  the  classic  waterfall,  but  here  there  were 
also  some  unexpected  instabilities  where  the  attachment  line  of  the  free 
surface  of  the  liquid  with  the  plate  developed  corrugations. 


INSTABILITY  AND  CHAOTIC  BEHAVIOUR  IN  A  FREE-SURFACE  FLOW 
W.  G.  Pritchard 

1 .  INTRODUCTION 

A  description  is  given  of  experimental  observations  of  a  free-surface 
flow  that  was  highly  unstable,  exhibiting  a  wide  range  of  different  motions, 
over  seven  hundred,  some  of  which  were  steady,  some  nearly  periodic  in  time 
and  some  of  which  had  a  less  regular,  almost  chaotic,  behaviour  in  time.  The 
particular  situation  to  be  studied  resembled  that  of  a  'waterfall1,  where 
fluid  flowed  down  the  bed  of  a  uniform  channel  and  poured  over  the  end  of  a 
plate  into  a  reservoir.  It  was  the  nature  of  the  flow  in  the  region  near  the 
end  of  the  plate  that  was  of  special  interest  here,  with  the  falling  liquid 
having  the  propensity  to  develop  many  different  flow  structures.  Because  the 
Reynolds  numbers  for  these  experiments  were  not  large,  ranging  between  0.05 
and  37,  it  would  appear  that  the  non-uniqueness  was  controlled  mainly  through 
the  nonlinearities  at  the  free  surface. 

Although  a  specification  of  the  flow  conditions  depended  on  at  least 
seven  dimensionless  parameter  groupings,  only  one  quantity,  namely  the  flux, 
was  systematically  varied  in  these  experiments.  At  very  small  values  of  the 
flux  the  flow  consisted  basically  of  a  number  of  sites  from  which  droplets  of 
liquid  fell  into  the  reservoir.  As  the  flux  was  increased  continuous  streams 
became  possible  and  at  larger  values  of  the  flux  the  liquid  had  a  tendency  to 
fall  in  the  form  of  sheets  rather  than  in  the  form  of  streams.  Various 
combinations  of  drops  and  streams  or  streams  and  sheet-flows  were  possible, 
leading  to  the  large  multiplicity  referred  to  above.  Some  of  the  observed 
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flow  structures,  such  as  the  (2T  +  C)  and  the  (2T  +  20  flows  described 


below,  had  a  rather  complicated  dependence  on  time.  In  these  cases  the 
temporal  dependence  was  not  reminiscent  of  a  nearly-periodic  motion  but  was  of 
a  more  chaotic  nature,  with  the  flow  differing  markedly  from  'cycle'  to 
' cycle' . 

The  paper  has  a  second  purpose,  concerning  the  position  at  which  the  free 
surface  of  the  liquid  attached  to  the  plate.  There  have  been  many  attempts  in 
recent  years  to  simulate,  by  numerical  computation,  certain  free-surface  flows 
that  arise  in  so-called  die-well  problems  and  in  the  coating  of  solid  surfaces 
with  liquid  films.  It  is  common  in  such  calculations  (e.g.  see  Nickell, 

Tanner  and  Caswell  1974,  Silliman  and  Scriven  1978,  Omodei  1980,  Cuvelier 
1981,  Saito  and  Scriven  1981,  Dutta  and  Ryan  1982)  to  specify  a  priori  the 
attachment  point  of  the  free  boundary  to  the  rigid  surface;  usually  a  sharp 
corner  of  the  rigid  surface  is  chosen  for  the  attachment  point.  However,  Jean 
and  Pritchard  (1980)  questioned  the  validity  of  this  assumption  and 
demonstrated  in  laboratory  experiments  that  it  does  not  hold  in  general.  It 
is  the  thesis  of  Jean  and  Pritchard  that,  without  the  means  of  pinning  the 
free-surface  attachment  at  some  particular  place  through  the  application  of 
'external'  point  forces,  the  attachment  position  of  the  free  surface  needs  to 
be  determined  as  part  of  the  overall  solution  of  the  flow  problem.  Thus, 
assuming  the  free  surface  attaches  to  the  solid  at  a  sufficiently  smooth  part 
of  the  boundary,  both  Jean  (1980)  and  Solonnikov  (1980)  were  able  to  demon¬ 
strate  the  existence  of  a  unique  solution,  in  two  dimensions,  to  a  mathemati¬ 
cal  problem  representing  a  certain  class  of  flows  of  the  above  kind.  An 
important  feature  of  these  theories  is  that  there  should  be  a  surface  tension 
acting  at  the  free  boundary  of  the  liquid  phase  for,  in  the  absence  of  surface 
tension,  the  basic  estimates  deriving  from  the  theory  for  elliptic  equations 


were  not  strong  enough  to  allow  an  existence  theory  to  be  established.  Some 
numerical  calculations  have  recently  been  reported  by  Tanner,  Lam  and  Bush 
(1984)  in  which  the  attachment  point  was  not  specified  a  priori. 

In  the  present  paper  empirical  evidence  will  be  given  to  support  the  case 
argued  by  Jean  and  Pritchard  (1980)  that  the  location  at  which  a  free  surface 
contacts  a  rigid  boundary  is,  in  gneeral,  determined  by  the  overall  flow 
conditions.  In  particular,  examples  will  be  given  showing  that,  for  given 
flow  parameters,  more  than  one  stable,  steady  flow  was  possible,  and  that 
these  flows  were  distinguishable  by  the  locations  of  the  attachment  points  of 
the  free  boundaries  to  the  rigid  surfaces.  The  attachment  points  were  not,  in 
general,  seen  to  coincide  with  a  corner  of  the  rigid  boundary.  So,  it  would 
appear  that  usually  it  is  the  attachment  point  which  adjusts  to  reflect  the 
global  conditions  rather  than  the  downstream  conditions  adjusting  in  response 
to  a  preferred  location  of  the  contact  point. 


2 .  EXPERIMENTAL  APPARATUS 


General  layout 

A  schematic  layout  of  the  apparatus  is  given  in  figure  1.  The  central 
part  of  the  apparatus  was  an  open  channel  of  uniform  section,  of  width  127  nun 
and  depth  74.5  mm,  sloping  at  an  angle  0.0735  rad  to  the  horizontal;  it  was 
1.20  m  long.  The  bed  of  the  channel  was  made  from  sheet  duralumin  plate,  of 
thickness  6.4  mm,  resting  on  an  alluminum  C-section  (not  indicated  in  the 
diagram).  At  the  upstream  end  of  the  channel  was  an  inlet  reservoir  from 
which  the  liquid  (an  oil)  poured  into  the  channel,  and  at  the  downstream  end 
was  a  receiving  reservoir.  The  side  walls  of  the  channel  and  of  the 
reservoirs  were  10  mm-thick  sheet  acrylic  (’Perspex'). 

The  return  from  the  reservoir  downstream  to  the  one  upstream  was 
basically  a  25  mm-bore  pipe.  A  T-joint  into  this  pipe,  just  below  the  exit 
from  the  downstream  reservoir,  led  through  a  stopcock  to  an  auxiliary 
reservoir.  This  reservoir  was  seated  on  a  platform  that  could  easily  be 
raised  or  lowered,  conveniently  allowing  liquid  either  to  be  added  to  or  taken 
from  that  which  was  circulating  round  the  main  system.  The  flow  was  forced  by 
a  centrifugal  pump,  located  as  shown  in  the  sketch.  The  pump  was  driven  by  a 
(commercially  produced)  d.c.  motor-tachogenerator  with  a  feedback  speed- 
control  system,  the  speed  of  the  motor  nominally  being  controlled  to  within  1 
part  in  104.  After  leaving  the  pump  the  oil  passed  through  a  heat  exchanger 
and  a  filter  before  returning  to  the  upstream  tank.  The  heat-exhange  unit 
consisted  of  an  array  of  copper  tubes  through  which  the  oil  passed  and  around 
which  cool  water  flowed.  By  carefully  adjusting  the  flow  rate  of  the  water 
past  the  array  of  tubes  a  nicely  controlled,  stable  operating  temperature 
could  be  achieved  for  the  circulating  oil.  A  perforated  plate,  which  was 
located  just  above  the  inlet  to  the  upstream  reservoir,  helped  even  out  the 
flow  in  the  reservoir. 
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A  mercury  bulb  thermometer  was  located  in  the  upstream  reservoir,  as 
shown  in  the  sketch  (figure  1).  The  oil  temperature,  which  was  noted  at 
frequent  intervals,  was  measured  to  an  accuracy  of  approximately  ±  0.01  °C. 

The  temperature  of  the  laboratory  was  maintained  at  (21  ±  0.5)°C  and  the  oil 
temperature  was  usually  stabilized  at  about  21°C. 

The  main  unit  rested  on  feet  located  under  each  of  the  reservoirs. 
Adjustment  screws  in  the  feet  were  carefully  set  to  ensure  that  the  bed  of  the 

channel  was  nearly  horizontal  in  a  direction  normal  to  the  main  axis  of  the 

channel  (i.e.  in  a  direction  normal  to  the  side  walls).  In  particular,  at  the 
end  of  the  plate  that  protruded  into  the  downstream  reservoir,  the  bed  was 

horizontal  in  the  cross-channel  direction  to  within  0.02  mm.  The  slope  of  the 

bed  in  the  long-channel  direction  was  determined  by  direct  measurement 
relative  to  the  surface  of  a  stationary  liquid  in  the  channel. 

More  complete  details  of  many  of  the  experimental  techniques  are  to  be 
reported  in  a  forthcoming  paper  concerning  free-surface  flows  of  liquids  in 
channels  at  small  Reynolds  numbers,  for  which  purpose  the  present  apparatus 
was  originally  constructed. 


The  downstream  reservoir 

The  experiment  centered  around  the  free-surface  patterns  observed  in  the 
downstream  reservoir,  the  details  of  which  are  given  in  figure  2.  The 
reservoir  was  basically  a  Perspex  box  of  rectangular  cross  section  into  which 
protruded  the  duralumin  plate  used  to  form  the  bed  of  the  main  channel.  The 
end  of  the  plate  was  carefully  milled  to  be  square  to  its  upper  and  lower 
surfaces,  and  edges.  The  plate  formed  an  'airtight'  fit  with  the  walls  of  the 
reservoir . 
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The  liquid  employed  in  these  experiments  was  a  commercial  vegetable  oil 
used  in  the  catering  industry.  This  oil  was  a  blended  mixture  of  a  variety  of 
vegetable  oils,  but  is  believed  predominantly  to  have  been  derived  from 
oilseed  rape.  One  minor  inconvenience  of  this  oil  was  that  its  mechanical 
properties,  in  particular  its  shear  viscosity,  showed  small  changes  with 
use.  This  kind  of  oil  has  been  used  for  about  five  years  in  our  laboratory 
and  over  this  period  it  has  been  found  that,  in  normal  use,  the  shear 
viscosity  of  the  oil  increased  by  roughly  0.04%  per  day.  Thus,  regular 
monitoring  of  the  viscosity  of  the  oil  was  needed.  This  was  usually  done  once 
a  week  during  periods  of  experimentation.  Some  tests  also  were  made  to 
ascertain  whether  or  not  the  oil  showed  significant  non-Newtonian  properties 
by  estimating  its  primary  normal-stress  difference  for  a  viscometric  flow. 
These  tests  were  carried  out  in  the  laboratory  of  Professor  A.  S.  Lodge  at  the 
University  of  Wisconsin  by  Mr.  H.  Vo.  The  liquid  indicated  negligible  normal- 
stress  effects,  even  at  shear  rates  several  orders  of  magnitude  in  excess  of 
those  arising  in  the  present  experiments.  Also,  over  the  same  range,  the 
shear  viscosity  was  found  to  be  independent  of  the  shear  rate. 

Near  the  usual  operating  temperatures  for  the  present  experiments  the 
(shear)  viscosity  of  the  oil  was  found  to  change  by  approximately  3.7%  per 
degree  c»  de ,  the  viscosity  decreasing  with  increasing  temperature. 

Therefore,  no  ily  was  it  important  to  monitor  carefully  the  temporal 
properties  c  r.e  oil,  hut  it  was  also  necessary  to  determine  the  operating 
temperature  for  each  experiment. 

The  viscosity  of  the  oil  was  monitored  by  measuring  the  time  for  a  given 
quantity  of  fluid  to  discharge  through  a  capillary  viscometer.  This  time  was 
determined  to  a  reproducible  accuracy  of  better  than  0.1%,  and  so  relative 
changes  of  the  viscosity  of  the  oil  could  be  tracked  very  closely.  The 


More  complete  details  of  these  flows  will  be  given  below  (in  particular, 
the  end  elevation  relating  to  figure  4(d)  reveals  additional  structure),  but 
the  photographs  of  figure  4  suffice  for  the  present  to  illustrate  the 
importance  of  allowing  the  attachment  point  to  adjust  itself  according  to  the 
global  flow  properties. 

' 

4.2.  The  overall  picture 

There  was  such  a  large  variety  of  possible  flows  in  this  experiment  that 
a  description  of  the  overall  picture  is  quite  complex.  For  this  reason  I 
shall  introduce  special  terminology,  as  the  description  unfolds,  to  be  able  to 
refer  easily  to  the  various  flow  structures. 

Five-site  flows 

At  the  smallest  flow  rates  that  could  be  maintained  in  these  experiments, 
the  flow  over  the  end  of  the  plate  settled  into  a  form  indicated  by  the  photo¬ 
graphs  of  figure  5.  The  oil  flowed  round  the  end  of  the  plate  to  the  under¬ 
side,  where  a  row  of  droplets  formed.  Eventually  the  droplets  would  acquire 
too  much  liquid  to  remain  attached  to  the  plate  and  would  fall  into  the  main 
reservoir,  at  which  stage  a  new  droplet  would  develop  at  the  site  just  vacated 
by  the  old  droplet.  For  this  particular  branch  of  flows  the  droplets  formed 
at  five  sites.  We  shall  refer  to  such  a  flow  as  a  ( 5 D )  (or  (D  D  D  D  D))  flow, 
to  indicate  that,  at  each  of  the  'sites'  from  which  liquid  fell  from  the 
plate,  the  motion  consisted  of  the  periodic  formation  of  droplets,  where  the 
above  sequence  of  D's  corresponds  to  the  sites  as  viewed  from  left  to  right  in 
the  photographs. 
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Photographs  were  also  taken  of  a  side  elevation  of  the  experiment,  the 
illumination  for  which  was  the  same  as  for  the  end-elevation  photographs . 

Thus,  bright  regions  of  the  photographs  indicate  that  light  rays  have  passed 
through  the  channel  to  the  camera  lens,  whereas  the  dark  regions  indicate  that 
rays  from  the  source  had  either  been  reflected  or  refracted  away  from  the 
camera  at  a  surface  of  the  liquid  sheet-  The  complex  structures,  associated 
with  the  stream  running  down  the  wall  of  the  channel  and  with  the  flow  near 
the  reservoir  walls  will  be  evident  from  some  of  these  photographs. 

Before  describing  the  different  flows  that  have  been  observed  in  these 
experiments  we  shall  first  illustrate  the  importance  for  theoretical  models  of 
not  specifying  a  priori  the  attachment  point  of  the  free  surface. 

4.1.  The  attachment  of  the  free  surface  to  the  plate 

In  figures  4(a),  (b)  are  shown  side  elevations  of  two  flows  made  under 
essentially  the  same  operating  conditions.  These  are  representative  of  two 
branches  of  flows  in  which  a  sheet  of  liquid  spanned  the  width  of  the  channel 
and  was  attached  to  the  plate  on  its  underside.  The  flows  in  figures  4(a), 
(b),  at  R  ~  9.2,  can  be  distinguished  by  the  distance  of  the  attachment 
point  from  the  corner  of  the  plate.  For  the  flow  shown  in  figure  4(a)  this 
distance  was  approximately  1.13a  whereas,  for  the  flow  shown  in  figure  4(b), 
the  attachment  point  was  only  about  0.32a  from  the  corner. 

Another  example  of  the  same  kind  of  phenomenon  is  shown  in  figures  4(c), 
(d).  The  flow  of  figure  4(c)  was  from  the  same  branch  of  flows  as  the  one 
shown  in  figure  4(b),  but  with  R  ~  23.1.  In  this  case  the  attachment  point 
on  the  underside  of  the  plate  was  approximately  0.15a  from  the  corner.  The 
photograph  of  figure  4(d)  at  R  ~  23.3  was,  however,  representative  of  a 
branch  of  flows  in  which  the  attachment  point  of  the  free  surface  was  on  the 
end  face  of  the  plate,  roughly  in  the  middle  of  this  face. 
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4.  EXPERIMENTAL  RESULTS 

The  discussion  of  the  experimental  observations  will  be  centered  around 
photographs  showing  the  main  features  of  the  visible  flow  properties.  The 
photographs  were  taken  from  one  of  two  aspects.  The  first  was  an  end 
elevation,  looking  upstream  along  the  channel.  For  these  photographs  the 
camera  was  aimed  slightly  downwards,  thus  rendering  visible  the  free  surface 
of  the  oil  in  the  reservoir  (and  sometimes  reflections  in  that  surface  are 
evident  in  the  photographs).  A  length  scale  can  be  taken  approximately  from 
the  plate  thickness  a,  but  note  that  the  plate  was  photographed  through  a 
curved  liquid  surface  which  acted  as  a  lens.  The  significance  of  this  effect 
can  be  seen  by  the  apparent  distortion  in  some  photographs  of  the  plate  near 
the  walls  of  the  channel.  Alternatively,  a  length  scale  is  obtainable  from 
the  channel  width  b,  though  this  can  be  difficult  to  discern  accurately  in 
some  of  the  photographs.  The  field  was  illuminated  from  the  right  (for  all 
photographs)  and,  in  order  to  provide  adequate  illumination  to  the  entire 
field,  the  stream  on  the  right-hand  wall  was  sometimes  over-exposed  in  the 
photographs.  Near  the  top  of  the  photographs  showing  the  end  elevation,  just 
above  the  plate,  is  a  'horizontal'  line:  this  is  a  reflection  from  the  curved 
surface  of  the  liquid  pouring  over  the  end  of  the  plate.  The  blurred  images 
in  the  background  of  these  photographs  are  some  screws  (used  in  the 
construction  of  the  reservoir)  and  the  end  of  the  aluminum  C-section  used  to 
support  the  channel. 

In  the  experiments  there  was  a  distinct  stream,  or  rivulet,  on  the  side 
walls  of  the  reservoir,  one  edge  of  which  connected  with  the  edge  of  the 
meniscus  on  the  wall  of  the  channel,  the  other  edge  forming  part  of  the 
meniscus  line  on  the  wall  of  the  reservoir  under  the  plate. 
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ai  Y-j  and  y ^  were  held  constant.  The  depth  h  of  the  approaching  stream 
was,  however,  varied  and  so,  in  the  parameterization  (3.9),  the  only 
parameters  to  remain  fixed  were  a,  Y1 ,  y 2  and  G.  The  reason  h  was  chosen 
as  the  length  scale  was  that  it  seemed  the  scale  most  relevant  to  the  balances 
associated  with  R,  s  and  G.  In  describing  the  experiments  we  shall, 
therefore,  quote  only  R,  S,  h  and  v,  from  which  all  other  parameter  values 
can  be  deduced,  including  the  velocity  scale,  if  required. 
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f(x,t)  =  0  , 

where  f  is  a  suitable  smooth  function  constrained  by  the  fact  that  points  in 
the  free  surface  remain  there  for  all  time.  This  leads  to  the  kinematic 
condition  on  f: 

ft  +  «  *  Vf  =  o  for  &  e  rD  .  0.6) 

Let  n(x)  and  t(x)  be  mutually  orthogonal  unit  vectors  at  each  point 
x  e  T  ,  with  n  normal  to  the  surface  and  directed  out  of  SI,  Then  the 

r  ^ 

boundary  conditions  at  the  free  surface  arising  from  dynamical  considerations 
are : 
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The  length  scale  which  most  appropriately  characterizes  ft  would  appear 
to  be  h,  and  a  representative  velocity  scale  would  appear  to  be  U  :=  (Q/h), 
with  Q  specified  by  (2.2).  Using  these  scales  we  shall  henceforth  consider 
all  variables  to  be  in  dimensionless  form. 

Let  lj(x)  denote  the  fluid  velocity  at  each  position  x  e  ft  and  let 
p(x)  denote  the  pressure.  Then  the  dynamical  equations  governing  the  fluid 
motion  in  ft  are 

9u 

R(-yjr  +  (u*V}u)  =  -Vp  +  Au  -  ,  (3.2) 

and  V«u  =  0  ,  (3.3) 

where  t  denotes  time,  j  is  the  unit  vector  in  the  vertically  upwards 

2 

direction,  R  :=  Uh/v  (=  Q/v  here)  and  G  :=  gh  /VU  (=  3/sin  a  here).  In 
these  definitions  v  is  the  kinematic  viscosity  of  the  fluid  and  g  is  the 
gravity  constant.  The  quantity  R,  the  Reynolds  number,  characterizes  the 
relative  importance  of  inertial  and  viscous  stresses ,  and  G  is  the  gravity 
parameter,  typifying  the  relative  importance  of  the  body  forces  due  to  gravity 
and  the  viscous  stress  field.  Note  that  (R/G)  -  U//gh  is  the  parameter 
usually  referred  to  as  the  Froude  number.  (The  gravity  parameter  G  took  the 
value  40.9  in  the  experiments,  F  ranged  between  0.04  and  0.95  and  the 
Reynolds  number  between  0.05  and  37.) 

The  stress  tensor  0(£),  which  is  scaled  by  (pvu/h),  where  p  is  the 
fluid  density,  is  given,  in  rectangular  cartesian  coordinates,  by 

°i j  :=  "P6ij  +  (uj,i  +  ui,j>  •  (3'4) 

The  boundary  conditions  applying  to  (2.1)  are  as  follows. 

(i)  x)  =  g  for  x  e  Td  ,  (3.5) 

where  g  is  some  prescribed  velocity  field  (compatible  with  (3.3)). 

(ii)  Let  us  suppose  that  the  set  T  representing  the  free  surface  is 
specified  implicitly  by  the  equation 


The  flow  structures  described  below  were  apparently  very  sensitive  to  changes 
in  d. 

Let  the  depth  in  the  channel,  well  upstream  of  the  end  of  the  plate,  be 
h. 

Then,  in  summary,  it  would  appear  that  the  parameters 

a,  b,  d,  a,  h  (3.1 ) 

suffice  to  fix  the  geometric  features  of  the  present  experiment  together  with 
(i)  the  distances  of  the  end  walls  of  the  reservoir  from  then  end  of  the 
plate,  (ii)  the  height  of  the  end  of  the  plate  above  the  floor  of  the 
reservoir  and  (iii)  the  location  and  size  of  the  outlet  hole  in  the  bottom  of 
the  reservoir. 

All  these  quantities,  with  the  exception  of  h,  were  kept  fixed  in  the 
experiments.  Values  for  the  parameters  mentioned  in  (i),  (ii)  and  (iii)  can 
be  obtained  from  figure  2,  but  it  is  believed  they  did  not  play  a  prominent 
role  in  the  outcome  of  the  experiments.  The  values  for  the  other  parameters 
were:  a  =  6.4  mm,  b  =  127  mm,  d  =  36  mm  and  a  =  0.0735.  Note  that, 

although  the  actual  value  of  d  was  not  determined  very  accurately,  it  was, 
however,  held  fixed  relative  to  a  mark  on  the  wall  of  the  reservoir,  to  much 
greater  accuracy,  as  described  above. 

We  now  determine  the  parameters  that  arise  from  dynamical  considera¬ 
tions.  The  discussions  of  §2  suggest  that  the  fluid  motions  should  be  well 
described  by  the  incompressible  form  of  the  Navier-Stokes  equation.  Suppose 
the  flow  occupies  an  open  domain  D  c  and  let  3fl  denote  the  boundary  of 

£5.  Let  3S2  be  partitioned  into  disjoint  sets  T  and  T  on  which 

Dirichlet-type  conditions  and  free-surface  conditions  are  to  be  specified 
respectively.  Note  that  ii,  and,  in  the  present  case,  T  may  depend  on 

r 

time. 
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FIGURE  3. 
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3 .  A  PARAMETERIZATION 

The  aim  here  is  to  provide  a  maximal  list  of  parameters  that  are  needed 
to  specify  the  flow  conditions.  Consider  a  flow  in  a  channel  of  width  b, 
the  bed  of  which  slopes  at  angle  a  to  the  horizontal.  We  wish  to  describe 
the  flow  field  in  the  region  near  the  end  of  a  plate  of  thickness  a,  which 
protrudes  into  the  downstream  reservoir.  There  the  fluid  pours  off  the  plate 
into  the  reservoir,  as  depicted  schematically  in  figure  3. 

In  some  of  the  experiments  to  be  reported  there  was  an  unbroken  sheet  of 
liquid  spanning  the  reservoir,  so  that  the  possibility  existed  of  sustaining  a 
differential  in  pressure  across  the  liquid  sheet.  We  shall,  however,  avoid 
the  complicating  factor  of  this  parameter  in  the  present  experiments  by 
ensuring  that  atmospheric  pressure  is  maintained  on  either  side  of  the  sheet 
( cf .  figure  3  )  . 

Suppose  that  the  lower  corner  of  the  end  of  the  plate  is  at  a  height  d 
above  the  tree-surface  level  in  the  reservoir  well  away  from  the  zone  in  which 
the  liquid  falls  into  the  reservoir.  That  d  could  be  defined  in  this  manner 
in  the  present  experiments  means  that  the  end  walls  of  the  reservoir  were 
effectively  at  'large  distances'  from  the  end  of  the  plate.  In  all  the 
experiments  below  d  was  kept  fixed,  in  the  following  sense.  Because  of  the 
fixed  volume  of  liquid  circulating  in  the  system  d  changed,  in  response  to 
the  changed  depth  of  liquid  flowing  down  the  channel,  whenever  the  flow  rate 
was  changed.  Compensation  for  the  consequent  changes  in  d  were  made  by 
adjustments  from  the  auxiliary  reservoir,  a  process  that  was  easy  to  carry  out 
in  the  present  experiments.  Usually  d  was  set  to  within  about  0.01  mm  of  a 
fixed  value  (by  bringing  the  free-surface  level  to  a  mark  on  the  wall  of  the 
reservoir)  and  usually  changes  to  the  flow  rate  were  small  enough  that  the 
adjustments  needed  to  keep  d  at  its  prescribed  value  were  less  than  0.5  mm. 
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error  in  the  specification  of  Q  (cf.  2.2).  The  temperature  of  the  oil  was 
determined  to  within  approximately  ±  0.02°C,  corresponding  to  an 
uncertaintity  of  ±  0.2%  in  the  nominal  value  of  v. 

It  will  be  important  later  to  have  some  idea  of  how  nearly  constant  Q 

remained  throughout  any  given  experiment.  As  indicated  above,  the  speed  of 

the  pump  used  to  drive  the  flow  was  controlled  by  a  feedback  system  which,  on 

the  basis  of  previous  experience  in  the  laboratory,  should  have  provided  a 

very  tight  control  on  the  flux.  Also,  to  check  that  the  depth  measurements 

did,  indeed,  give  an  accurate  representation  of  the  flux,  some  direct 

measurements  were  made  of  surface  speeds  in  the  channel.  These  measurements 

showed  less  than  ±  0.5%  variation  in  surface  speed  at  each  of  a  number  of 

pump  settings,  suggesting  that  the  flux  was,  in  fact,  held  very  nearly 

constant.  On  the  other  hand,  depth  measurements  made  under  ostensibly  steady 

flow  conditions  showed  occasional  deviations  of  up  to  0.02  mm.  Measured  depth 

variations  of  this  magnitude  were  significant  and  correspond,  in  a  depth  of 

say  3  mm,  to  a  2%  change  in  Q.  I  do  not  know  how  such  variations  came  about 

as  small  effects  of  this  kind  could  arise  from  many  possible  sources:  for 

example,  the  sensitivity  of  the  system  to  small  effects  is  illustrated  by  the 

,  1  o 

feature  that  a  change  of  temperature  of  —  c  would  have  brought  about  a  depth 
change  of  approximately  0.3%  (or  approximately  0.01mm  in  a  depth  of  3mm),  not 
even  allowing  for  any  changes  in  the  pump  characteristics;  also,  there  is  a 
potential  long-wave  instability  of  the  free  surface  when  the  Reynolds  number 
exceeds  a  critical  value  Rc  =  -g-  cot  a  (e.g.  see  Yih  1969,  p.  502)  but,  as 
many  of  the  measurements  were  made  at  Reynolds  numbers  well  below  the  value  of 
11.3  for  Rc  obtaining  in  the  present  experiment,  this  effect  would  not 
appear  to  provide  the  explanation  for  the  observed  depth  variations. 
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for  steady  flow  down  the  plane  is 

2 

gh  sin  a  (  rVl2 I 

u  =  “  Jv  i2  (p  ~  ®  J  '  (2.D 

where  v  is  the  kinematic  viscosity  of  the  fluid  and  g  is  the  gravity 
constant.  Corresponding  to  this  flow,  the  volume  flux,  Q,  per  unit  width  is 

Q  =  (gh3  sin  a)/3v  ,  (2.2) 

and  the  velocity  at  the  free  surface  is  equal  to  (1.5  Q/h).  Thus,  the  liquid 
depth  in  the  uniform  channel  provides  a  convenient  means  of  specifying  the 
basic  flow  parameters  for  the  present  experiment.  That  the  above  distribution 
(2.1)  did  indeed  give  a  good  description  of  the  flow  in  the  channel  was 
checked  by  direct  measurements  of  the  surface  speed  which  agreed,  over  a  range 
of  operating  conditions,  to  within  1.5%  of  that  predicted  from  (2.2)  using  the 
empirically  determined  values  of  h,  a  and  v. 

Measurements  of  the  depth  h  were  made  near  the  central  plane  about 
midway  along  the  channel  using  a  pointer  gauge.  (Except  for  small  regions 
near  the  walls  and  near  the  entrance  and  exit  of  the  channel,  the  depth  was 
found  to  be  very  nearly  uniform  throughout.)  This  gauge  consisted  of  a  fine 
sewing  needle  mounted  in  the  end  of  a  piece  of  1.2  mm-diameter  brass  rod  which 
could  be  raised  and  lowered  in  a  direction  normal  to  the  plane  of  the  channel 
bed.  A  dial  gauge  bearing  on  the  upper  end  of  the  brass  rod  was  used  to 
determine  the  height  of  the  needle  point  above  the  bed  of  the  channel.  The 
dial  gauge  was  graduated  in  units  of  0.01  mm  and  had  a  range  of  25  mm.  The 
position  of  the  free  surface  was  found  by  carefully  lowering  the  needle  until 
the  point  broke  the  surface,  an  event  that  was  easily  recognizable  because  of 
surface  tension  pulling  the  liquid  up  the  shaft  of  the  needle.  The  depth 
measurements  were  accurate  and  reproducible  to  within  approximately  ±  0.01 
mm.  Thus,  for  a  depth  of  1  mm,  corresponding  roughly  to  the  smallest  flow 
rate  used  here,  the  accuracy  of  the  determination  of  h  corresponds  to  a  3% 
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absolute  viscosity  was  determined  by  comparison  of  the  viscosity  of  the  oil 
with  that  of  distilled  water,  the  value  of  which  was  taken  from  international 
standards-  The  kinematic  viscosity  of  the  oil  at  the  working  temperatures  of 
the  experiment  was  approximately  0.8  Stokes  (or  80x10  ^  m^s-*),  which  is 
roughly  80  times  the  kinematic  viscosity  of  water  at  the  same  temperature. 

The  actual  viscosity  for  a  given  experiment  was  known  to  within  1%. 

The  surface  tension  of  the  oil  was  determined  in  two  ways.  The  force 
exerted  on  a  thin  plate  whose  edge  had  been  brought  down  to  touch  the  surface 
of  the  oil  was  measured,  and  this  quantity  was  used  to  estimate  the  suface 
tension.  (Measurements  were,  of  course,  made  with  plates  of  different  lengths 
to  allow  for  the  anomalies  at  the  ends  of  the  plate.  Also  the  contact  angle 
of  the  oil  with  the  plate  surface  was  observed  to  be  very  small,  if  not  zero.) 
The  second  method  consisted  of  determining  the  force  acting  on  a  circular  ring 
of  wire  drawn  out  through  the  liquid  surface. 

From  these  two  kinds  of  measurements  the  surface  tension  of  the  oil 
(determined  at  a  temperature  of  21°C)  was  estimated  to  be  33.0  dyn/cm  ±  3%. 

The  density  of  the  oil  was  measured  using  a  relative  density  bottle  and 
was  estimated  to  be  0.922  gm  cm-3  at  a  temperature  of  21°C. 

Specification  of  the  flux 

The  basic  physical  quantity  being  varied  in  this  experiment  was  the 
volume  flux  through  the  system.  A  convenient  way  of  monitoring  this  quantity 
is  to  measure  the  depth  of  the  liquid  flowing  down  the  channel. 

Consider  a  viscous  fluid  flowing,  under  the  action  of  gravity,  down  a 
plane  surface  inclined  at  an  angle  a  to  the  horizontal.  Suppose  that  the 
thickness  of  the  sheet  of  liquid  is  h  and  that  y  measures  the  normal 
distance  above  the  plane.  Then  the  Poiseuille  velocity  distribution  u(y) 


The  shutter  speed  for  these  photgraphs  was  usually  (1/120)  s,  during 
which  time  the  droplets  could  fall  a  significant  distance,  accounting  for 
their  blurred  outline  :..i  some  of  the  photographs.  The  side  elevation  in 
figure  5(a)  provides  another  view  of  the  structure  of  this  branch  of  flows, 
and  it  can  be  seen  that  the  droplets  formed  at  a  distance  of  roughly  0.7a  from 
the  lower  corner  of  the  plate.  The  photograph  captured  a  droplet  just  at  the 
beginning  of  its  fall;  the  oil  level  in  the  reservoir  can  be  seen  near  the 
bottom  of  the  photograph. 

The  end-elevation  photographs  shown  in  figures  5(b)  -  (d)  are  representa¬ 
tive  of  this  branch  of  flows.  The  photograph  in  figure  5(b),  for  which 
R  ~  0.06,  was  taken  at  one  of  the  smaller  flow  rates  obtainable  in  the  current 
apparatus.  Although  it  may  appear  from  the  photograph  not  to  be  the  case, 
this  was  in  fact  a  (5D)  flow:  the  second  site  (from  the  left)  shows  the 
residual  string  of  liquid  after  the  droplet  had  fallen:  this  string  would 
have  'nipped'  itself  at  the  neck,  just  under  the  plate,  shortly  after  the 
photograph  was  taken.  A  circular  pattern  of  waves,  emanating  from  the  entry 
of  the  droplet,  can  be  seen  in  the  liquid  surface  in  the  reservoir.  (The 
small  air  bubbles  visible  in  the  reservoir  were  formed  by  the  falling 
droplets.)  The  { 5 D )  branch  of  flows  depicted  in  figures  5(b)  -  (d)  lost 
stability,  at  a  flow  rate  equivalent  to  R~  0.19,  to  the  ((5D)C)  flow  shown 
in  figure  5(e)  -  a  flow  consisting  of  5  drip  sites  and  one  site  on  the  extreme 
right-hand  side  that  was  a  continuous  (C)  stream.  (Note  the  visual  difference 
between  the  continuous  stream  in  figure  5(c)  and  the  dripping  flow  at  site  2 
in  figure  (5(b). I 

For  the  flows  at  the  'D'  sites,  the  droplets  formed  at  very  nearly 
regular  intervals,  the  period  of  which  decreased  as  the  flow  rate  increased. 
For  example,  in  one  sequence  of  experiments  the  formation  period  of  the  drop¬ 
lets  was  approximately  1.0  s  at  R  =  0.072,  decreasing  to  0.81  s  at 
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R  =  0.083 


to  0.68  s  at  R  =  0.95  and  to  0.52  s  at  R  =  0.144.  (In  that 


particular  experiment  the  (5D)  flow  lost  stability  to  a  ((5D)CC)  flow,  i.e.  a 
seven-site  flow,  at  R  ~  0.173.)  Although  the  droplet  formation  was  very 
nearly  a  periodic  process  it  was  not,  in  fact,  exactly  periodic  in  time,  a 
feature  that  was  apparent  because  the  phase  of  formation  of  the  drops  relative 
to  one  another  could  be  seen  to  change  gradually  with  time.  This  raises  a 
general  theoretical  question  of  whether,  in  the  absence  of  the  inevitable 
imperfections  in  a  laboratory  experiment,  the  flow  would  have  been  strictly 
periodic  in  time.  Certainly  it  was  my  impression  from  observing  the  motions 
that  a  strictly  period  flow  was  being  'spoiled'  by  imperfections,  resulting  in 
a  motion  that  was  just  nearly-periodic  in  time. 

The  (5d)  flow  described  above  was  the  most  stable  flow  observed  at  very 
small  flow  rates.  Indeed,  at  the  very  smallest  flow  rates  used  it  was  the 
only  flow  observed  in  these  experiments,  though  whether  or  not  this  branch  of 
flows  would  have  extended  down  to  the  limit  of  zero  flow  rate,  or  whether  it 
formed  a  unique  branch  at  small  flow  rates  is,  of  course,  not  known.  The 
particular  flow  to  which  the  (5D)  flow  lost  stability,  as  the  flow  rate  was 
increased,  was  quite  sensitive  to  the  'local'  conditions.  Thus,  in  three 
different  experiments  (two  of  which  are  mentioned  above)  the  (5D)  flow  lost 
stability  as  follows:  (i)  at  R  ~  0.17  to  a  ((5D)CC)  flow;  (ii)  at  R  ~  0.18 
to  a  ( ( 5 D ) C )  flow;  and  (iii)  at  R  ~  0.13  to  a  (6D)  flow.  In  a  separate 
experiment  a  (5D)  flow  was  established  when  a  six-site  flow  lost  stability 
at  R  ~  0.07. 


and 


Six-site  flows 

At  flow  rates  equivalent  to  Reynolds  numbers  between  approximately  0.068 
0.250,  branches  of  motions  characterized  by  six  sites  of  attachment  at  the 
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plate  were  possible.  An  example  of  one  such  motion  was  given  in  figure  5(e) 
where  a  (5D)  flow  lost  stability  to  a  ((5D)C)  flow,  the  latter  consisting  of 
five  drip  sites  and  one  site  that  was  a  continuous,  unbroken  stream.  I  was 
not  able  to  induce  the  formation  of  such  unbroken  streams  with  five-site 
flows,  but  they  could  easily  be  sustained  with  the  six-site  flows.  For 
example,  in  one  experiment  a  (5D)  flow  lost  stability  at  R  ~  0.13  to  the 
(6D)  flow  shown  in  figure  6(a).  Then  at  R  ^  0.162  a  variety  of  six-site 
flows,  typified  by  figures  6(b)  -  (h)  were  induced.  Thus,  at  these  operating 
conditions,  flows  comprising  one  unbroken  stream  and  five  drip  sites  could  be 
established  with  the  unbroken  stream  at  any  of  the  six  sites.  Examples  are 
shown  in  figures  5(b)  -  (d)  where  the  'C'  stream  was  at  sites  number  3,  5  and 
1  respectively.  Similarly,  flows  with  two  unbroken  streams  were  stable 
possibilities  at  these  conditions,  as  shown  by  figures  6(e)  -  (h).  (Note 
that,  in  figure  6(g)  sites  number  5  and  6  were,  in  fact,  drip  sites.)  These 
examples  show,  respectively,  (CC(4D)),  (CDC(3D)),  (CDDCDD)  and  ( DCDCDD)  flows, 
but  all  combinations  with  two  unbroken  streams  seemed  to  be  viable  states. 
Atttempts  to  establish  a  flow  under  these  conditions  with  three  unbroken 
streams  were,  however,  unsuccessful,  leading  in  all  cases  to  a  loss  of 
stability  of  the  six-site  flows  to  a  seven-site  flow.  In  the  experiments  it 
was  very  easy  to  transform  a  dripping  site  into  a  continuous  stream,  and  vice- 
versa,  by  simply  disturbing  the  flow  in  the  zone  on  the  upper  side  of  the 
plate  immediately  above  the  site  in  question.  This  was  usually  done  by 
dipping  a  glass  rod  of  about  4  mm  diameter  into  the  stream. 

Another  sequence  of  observations  concerning  six-site  flows  is  shown  in 
figure  7.  This  sequence  began  with  the  formation  of  the  ((5D)C)  flow  shown  in 
figure  5(e)  arising,  at  R  =  0.187,  from  an  exchange  of  stabilities  with  a 
(5D)  flow.  The  flow  rate  was  then  reduced  to  a  value  equivalent  to  R  =  0.150, 
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FIGURE  6.  Examples  of  nix-site  flowB.  (a)  (6D)  flow  at  p  ■=  O.lJh,  S  -  65 . 2 , 

h  •  1.51  mm,  v  ■=  0.700  St.  The  conditions  for  lh)-(h)  were:  R  '  0.162, 
S  ’  58.1,  h  -  1.60  run,  v  =  0.781  st .  (b)  (DDCMDI)  flow;  (c)  (MMCD) 

flow;  (d)  (CI5DM  flow;  (e)  (<"('(41)1)  flow;  (f)  (CU'(JD))  flow;  <  q  ) 
(CDICl)b)  flow;  (h)  (nCtK'DU)  flow. 


under  which  conditions  a  (6D)  pattern  (see  figure  7(b))  was  stable,  as  were 
patterns  comprising  one  'C'  stream,  two  'C*  streams  (see  figure  7(c))  and 
three  'C'  streams  (see  figure  7(d)).  The  only  pattern  with  three  'C'  streams 
I  established  here  was  the  one  shown  in  figure  7(d),  and  this  seemed  to  have  a 
rather  small  domain  of  stability.  I  was  unable  to  establish  a  six-site  flow 
with  four  *C'  streams.  Some  of  the  stability  properties  of  these  flows  may, 
however,  have  been  associated  with  the  feature  that  here  the  flux  Q  (and 
hence  R)  seemed,  from  the  depth  measurements  to  vary  by  bout  1  or  2%  over  a 
two-  to  three-minute  interval.  In  the  process  of  adjusting  the  flow  rate  from 
a  value  equivalent  to  R  =  0.187  to  R  =  0.142  a  stable  (DDCCDD)  flow 
changed  to  the  (6D)  flow  shown  in  figure  7(a).  The  (6D)  branch  of  flows  held 
stability  down  to  a  flow  rate  equivalent  to  R  ~  0.067,  at  which  point  a  (5D) 
flow  established  itself.  The  flow  rate  was  then  readjusted  such  that  R  =  0.162 
and  a  six-site  flow  was  reestablished.  At  these  conditions  I  induced  the  (6D) 
flow  and  the  'one-C'  and  'two-C'  flows  shown  in  figures  7(e),  (f),  (g),  but  I 
was  unable  to  induce  a  stable  pattern  with  more  than  two  continuous  streams. 
Among  the  possible  six-site  flows  at  R  =  0.202,  only  a  (6D)  pattern  (see 
figure  7(h))  appeared  to  be  a  stable  possibility,  with  attempts  to  induce  any 
'C'  streams  leading  to  a  loss  of  stability  of  the  six-site  flow.  The  (6D) 
pattern  naturally  lost  stability  at  R  ~  0.251  to  the  ( ( 7C ) D )  flow  shown  in 
figure  7 ( i  )  . 

All  the  six-site  flows  observed  were  nearly  periodic  in  time,  in  the  same 
sense  as  described  for  the  five-site  flows.  There  appeared  to  be  at  least  23 
different  six-site  flows  i>ossible  in  this  experiment. 
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Seven-site  flows 


A  class  of  flows  emanating  from  seven  sites  on  the  underside  of  the  plate 
were  observed  at  flow  rates  equivalent  to  R  ranging  between  approximately 
0.12  and  1.8.  Examples  of  some  of  the  flows  possible  under  these  conditions 
are  given  in  figures  8(a)  -  (p),  all  of  which  were  obtained  under  the  same 
operating  conditions;  many  more  such  flows  were  also  stable  under  the  same 
conditions.  Shown  in  figure  8(a)  is  a  flow  in  which  all  seven  sites  were  drip 
sites  (sites  3  and  4  were,  in  fact,  drip  sites),  whereas  (b)  and  (c)  give 
examples  in  which  one  of  the  sites  was  a  continuous  stream  -  namely  site  2  in 
figure  8(b)  and  site  7  in  (c)  -  but  the  C  stream  could  be  established  at  any 
of  the  seven  sites.  Examples  of  flows  involving  two  C  sites  are  shown  in 
figures  8(d)  -  (g).  Again,  all  combinations  involving  five  D  sites  and  two  C 
sites  seemed  to  be  possible.  It  also  seemed  possible  to  generate  any 
combination  involving  three  C  sites  and  four  D  sites  (some  examples  of  which 
are  shown  in  figures  8(h)  -  (k)),  four  C  sites  and  three  D  sites  (see  figures 
8(1),  (m),  ( n ) )  or  those  with  five  C  sites  and  two  D  sites  (e.g.  see  figure 
8(o)).  I  am  not  sure  if  flows  with  one  drip  site  and  six  C  sites  were  stable 
possibilities  under  these  operating  conditions  but  I  would  suspect  this  to 
have  been  the  case  for,  as  shown  in  figure  8(p),  a  steady  flow  with  seven  C 
sites  was  a  stable  possibility.  The  latter  is  the  first  steady  motion  so  far 
encountered  in  this  experiment,  all  the  other  flows  being  nearly  periodic  in 
time,  as  discussed  above.  Although  a  flow  with  one  drip  site  had  not 
specifically  been  induced  under  the  conditions  obtaining  for  figure  8,  such  a 
flow  was  observed  under  the  conditions  R  =  0.200,  S  =  52.7,  h  =  1.68  mm,  v  = 
0.755  St.;  I  was,  however,  unable  to  induce  the  steady  (7c)  flow  under  those 
conditions.  The  flows  involving  either  six  or  seven  C  streams  seemed  to  have 
quite  small  domains  of  stability,  though  no  specific  attempt  was  made  to 
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delineate  these  particular  domains  in  the  experiments.  But  estimates  were 
obtained  for  the  stability  limits  of  seven-site  flows  in  general.  So,  for 
example,  on  carefully  reducing  the  flow  rate  it  was  found  that  a  (D(5C)D)  flow 
lost  stability  to  a  (6D)  flow  at  R  =  0.128,  S  =  65.2,  h  =  1.51  mm,  v  =  0.804 
St.,  and  other  seven-site  flows  seemed  not  to  be  stable  possibilities  at 
these  conditions.  Near  the  upper  limit  it  was  found  that  the  (D(5C)D)  flow 
lost  stability  to  an  eight-stream  flow  when  R  =  0.265,  S  =  40.3,  h  =  1.92  mm, 
v  =  0.803  St.,  under  which  conditions  I  was  unable  to  reestablish  any  of  the 
seven-site  flows  by  disturbing  the  flow  near  the  end  of  the  plate. 

As  the  flow  rate  was  further  increased  beyond  the  stability  range  of  the 
above  flows,  a  different  kind  of  unsteady  motion  involving  seven  sites  was 
observed.  This  particular  flow  arose  when  I  was  attempting  to  establish  a 
nine-site  flow  from  an  eight-site  flow  and,  as  shown  in  figure  9,  it  consisted 
of  a  D  site  together  with  six  C  sites.  A  distinguishing  feature  of  this 
motion  was  that,  as  well  as  being  unsteady  by  virture  of  the  drip-formation 
process,  the  locations  of  the  continuous  streams  continually  changed  with 
time,  as  indicated  by  the  differences  between  figures  9(a)  and  (b).  Note 
that,  in  the  previous  examples  involving  both  D  and  C  sites,  a  small  adjust¬ 
ment  of  the  locations  of  the  d  streams  near  a  D  site  would  take  place  as  the 
droplet  grew  and  then  broke  off.  In  the  present  case,  however,  the  roughly 
periodic,  large-scale  excursions  of  the  C  streams  seemed  to  take  place  on  a 
time  scale  unrelated  to,  and  much  larger  than,  the  period  for  the  formation  of 
the  droplets.  The  impression  given  of  this  flow  was  that  a  'conventional' 
(D(6C) ]  flow  was  nearly  stable,  but  not  quite,  and  so  there  was  a  readjustment 
of  the  streams  to  allow  the  formation  of  a  new  stream;  however,  this  new 
stream  did  not  eventuate  and  the  flow  tried  to  equilibrate  itself,  through  the 
influence  of  surface  tension  forces,  so  that  the  spacing  of  the  streams  was 


-32- 


i  lLOkt  9.  An  u;isi  eady  (D(GC))  llow  ul  I-  0.294,  .*»  =  J7.2,  h  -  2.U0  nun,  v  0.009  St 


(a)  and  (b)  -ire  photographs  of  the  same  flow  at  two  different  times. 


nearly  uniform.  And  so  the  process  repeated  itself.  This  motion  was  observed 
for  over  16  minutes  and,  as  it  appeared  to  be  quite  stable,  the  operating 
conditions  were  then  changed.  Other  examples  of  this  kind  of  phenomena  will 
be  described  below. 

Another  branch  of  steady  (7C)  flows,  apparently  not  connected  to  the 
branch  corresponding  to  the  flow  shown  in  figure  8(p),  was  found  at  flow  rates 
corresponding  to  values  of  R  ranging  between  R  ~  0.79  and  R  ~  1.80. 
Examples  of  the  flows  observed  in  this  branch  are  shown  in  figures  10(a)  -  (c) , 
the  flow  losing  stability  near  R  =  1.80  to  the  one  shown  in  figure  10(f),  to 
be  referred  to  as  a  (T)  flow.  The  way  this  branch  lost  stability  at  the  lower 
end  of  the  range  is  discussed  below. 

Remark :  Although  the  two  branches  of  steady  (7C)  flows  shown  above  were 

not  connected  in  the  particular  region  of  parameter  space  investigated  here, 
these  branches  may,  of  course,  have  been  connected  via  some  other  path  through 
parameter  space . 

Some  side  views  of  the  present  branch  of  steady  (7C)  flows  are  given  in 
figure  11.  It  is  seen  that  the  attachment  line  of  the  free  surface  on  the 
underside  of  the  plate  was  approximately  0.9a  from  the  lower  corner  of  the 
plate.  The  place  at  which  the  jets  penetrated  the  free  surface  of  the 
reservoir  can  be  located  from  the  'dimple'  visible  at  the  free  surface  of  the 
reservoir . 

The  way  in  which  the  branch  of  steady  (7C)  flows  depicted  in  figures  10 
and  11  lost  stability,  as  the  flow  rate  was  reduced,  was  not  decisively  clear 
in  the  experiments.  It  is,  however,  my  view  that  stability  was  exchanged  for 
an  unsteady,  chaotic  type  of  motion,  comprising  basically  seven  continuous 
streams,  examples  of  which  are  given  in  figure  12.  The  reason  for  the 
difficulty  in  interpreting  the  experimental  observations  lies  in  the  possible 
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referred  to  as  (T)  flows. 


Further  examples  of  these  flows,  together  with  some 


side  views,  are  given  in  figure  18.  The  flow,  in  the  region  below  the  plate, 
consisted  of  two  nearly-cylindrical  jets  of  fluid  emanating  from  the  underside 
of  the  plate  and  connected  by  a  thin  sheet  of  liquid.  These  jets,  or  arms, 
were  pulled  towards  each  other  by  the  tension  in  the  liquid  sheet  and  it  can 
be  seen  from  the  figure  how  the  arms  became  more  nearly  vertical  as  they 
carried  a  larger  volume  of  liquid.  The  side  elevations  (figures  18  (d),  (f), 
(i),  (k))  show  how  the  attachment  line  of  the  free  surface  on  the  plate  was  on 
the  underside  of  the  plate,  well  behind  the  lower  corner  of  the  plate. 

The  branch  of  steady  (T)  flows  lost  stability,  with  increasing  flow  rate, 
at  R  =  4.65,  S  =  6.09,  h  =  4.94  mm,  v  =  0.790  St.  to  the  steady  flow  shown  in 
figure  18  (1),  to  be  referred  to  as  a  (QQ)  flow,  consisting  of  a  sheet  of 
liquid  spanning  the  channel  except  for  two  small  quadrants  missing  near  the 
junction  of  the  free  surface  in  the  reservoir  and  the  walls.  Note  that  the 
(T)  flow  shown  in  side  elevation  in  figure  18  (k)  was  at  a  slightly  larger 
value  of  R  than  that  at  which  the  branch  lost  stability.  The  reason  for 
this  was  that  the  viscosity  of  the  liquid  differ  substantially  from  that  which 
obtained  for  figures  18  (j)  and  (1),  suggesting  how  dependent  the  stability 
limits  can  be  on  changes  in  the  various  parameters. 

As  the  flow  rate  was  decreased  the  branch  of  steady  (T)  flows  lost 
stability  to  an  unsteady  motion  at  approximate  parameter  settings  of  R  = 

1.54,  S  =  12.7,  h  =  3.42  mm,  v  =  0.791  st.  The  photograph  shown  in  figure  19 
(d),  which  was  taken  at  these  conditions,  shows  the  pattern  at  a  reasonably 
quiescent,  and  nearly  symmetric,  phase  of  its  unsteady  cycle.  The 
unsteadiness  here  took  the  form  of  mild  oscillations  (in  the  plane  of  the 
sheet)  of  the  left-hand  arm  of  the  sheet;  these  oscillations  built  up  slowly 
in  amplitude  and  then  died  away.  At  a  slightly  reduced  flow  rate,  equivalent 
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Nine-site  flows 


[ 


Flows  emanating  from  nine  sites  on  the  underside  of  the  plate  were 
observed  at  flow  rates  corresponding  to  values  of  R  between  approximately 
0.18  and  1.4.  None  of  the  flows  in  this  class  arose  naturally  from  a  loss  of 
stability  of  another  flow.  They  were  all  established  by  making  suitable 
perturbations  to  the  stream  in  the  zone  above  the  end  of  the  plate.  As  with 
the  previous  examples,  combinations  of  drip  sites  and  unbroken  streams  were 
possible  though,  in  this  case,  it  was  slightly  more  delicate  than  hitherto  to 
establish  the  drip  sites.  In  particular,  I  was  able  to  induce  only  flows 
having  four  or  fewer  D  sites,  some  examples  of  which  are  given  in  figure  16. 
It  did,  however,  appear  as  though  all  combinations  of  such  flows  were  stable 
possibilities  and  it  was  my  impression  that,  given  more  patience,  I  might  have 
been  able  to  establish  some  flows  with  five  D  sites. 

I  was  unable  to  induce  drip  sites  with  this  class  of  flows  at  flow  rates 
above  those  corresponding  to  the  conditions  obtaining  for  figure  16  (for  which 
R  =  0.269),  and  it  seemed  as  though  no  nine-site  flows  were  stable  at  a  flow 
rate  below  that  equivalent  to  R  =  0.178. 

The  branch  of  steady  (9C)  flows,  some  examples  of  which  are  given  in 
figure  17,  was  stable  over  a  fairly  wide  range  of  flow  rates.  On  reducing  the 
flow  rate,  stability  of  the  branch  was  lost  at  R  =  0.178,  S  =  52.7,  h  =  1.68 
mm,  v  =  0.802  St.  to  a  steady  (8C)  flow.  The  steady  (9C)  branch  of  flows  also 
lost  stability,  on  increasing  the  flow  rate,  to  a  steady  (80  flow,  as  shown 
in  figure  17  ( f ) . 

T  flows 

Both  the  steady  ( 7C )  branch  of  flows  and  the  steady  (80  branch  of  flows 
lost  stability  (cf.  figures  10(f)  and  15(1))  to  a  branch  of  steady  flows  to  be 
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FIGURE  15  (ctd).  (h)  (80  flow:  R  *  1.22,  S  “  14.7,  h  “  3.18  on,  v  *  0.798  St.i  (1)  (80  flow:  R  “  1.53,  S  “  12.6, 
h  *  3.43  am,  V  -  0.797  St.i  (j)  (80  flow:  R  «  1.78,  S  “  11.5,  h  -  3.59  nun,  v  "  0.792  St.i  (k)  (6C) 
flow:  R  -  2.09,  S  -  10.4,  h  -  3.79  no,  v  «  0.790  St.r  (1)  The  branch  of  flows  lost  stability  to  this 


FIGURE  15.  Examples  from  the  steady  (BC)  branch  of  flows.  (a)  The  branch  of  flows  lost  stability  to  this  (6d)  flow: 

P  -  0.139,  F  «  63.5,  h  -  1.53  mm,  y  -  0.780  St .  j  <b)  (80  flow:  R  «  0.174,  S  -  54.6,  h  -  1.£5  «n,  v  -  0.7R9 

St. »  (c)  (80  flow,  pide  view:  R  -  0.260,  S  -  44.4,  h  •  1.83  mm,  y  -  0.7S4  St .  j  (d)  (RO  flow,  side  v1swt 

R  -«  0.498,  S  -  29.  1,  h  -  2.26  mm,  y  -  0.755  St-i  (•)  (BC)  flow:  R  -  0.559,  S  -  24.4,  h  -  2.47  mm,  v  -  0.806 

St.r  If)  <8C)  flow,  side  view:  R  •  0.741,  6  -  22.0,  h  •  2.60  mm,  v  •  0.756  St . »  (g)  (Hr)  flow:  R  •  0.863, 


v  =  0.806  St.,  all  combinations  of  flows  involving  zero  to  eight  drip  sites 
seemed  to  be  possible.  Then,  on  changing  the  operating  conditions  to  R  = 
0.282,  S  =  39.9,  h  =  1.93  mm,  v  =  0.784  St.,  I  was  unable  to  establish  eight- 
site  flows  with  more  than  four  drip  sites;  and,  under  the  conditions  R  = 
0.294,  S  =  37.2,  h  =  2.00  mm,  v  =  0.809  St.,  none  of  the  above  flows 
involving  a  drip  site  seemed  to  be  a  stable  possibility.  At  these  latter 
operating  conditions  there  was,  however,  one  unusual  flow  which  had  a  drip 
site  near  the  left-hand  wall  of  the  channel  and  seven  C  sites.  The 
difference  between  this  and  the  previous  'D  flows'  lay  in  the  rhythm  of  the 
dripping  process,  which  consisted  of  a  short  interval  between  two  drips, 
followed  by  a  longer  interval  until  the  formation  of  the  next  drip,  whereas 
normally  the  interval  between  successive  drips,  at  a  given  site,  was  very 
nearly  uniform.  This  atypical  D-flow  seemed  to  be  associated  with  asymmetries 
of  the  apparatus,  for  X  was  not  able  to  generate  a  flow  with  a  D  site  near 
the  right-hand  wall  of  the  channel  under  the  same  operating  conditions. 

Examples  spanning  the  range  of  flow  rates  for  which  the  branch  of  steady 
(8C)  flows  were  observed  are  given  in  figure  15.  On  reducing  the  flow  rate 
(cf.  figures  15  (g),  (e),  ( b ) )  the  branch  lost  stability  to  the  (6D)  flow 
•hown  in  figure  15  (a)  and,  on  increasing  the  flow  rate  (cf.  figures  15  (h)  - 
(k)),  the  (80  flow  eventually  lost  stability  to  the  (T)  flow  shown  in  figure 
15  (1).  Some  typical  side  views  of  the  flow  are  given  in  figures  15  (c),  (d), 
(f). 

All  eight-site  flows  observed  in  these  experiments  were  either  steady 
motions  (the  (80  branch)  or  nearly  periodic  in  time,  including  the  atypical 
(D(7C) )  flow. 
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FIGURE  14.  Ey.amjlea  of  flew-*.  (a)  f«D)  flow;  fS)  f(30)C(4D))  flow,  (c)  (OrKPrcn)  flew,  (d)  (DCUC'POCD) 

flow,  («)  (DCnur*CD)  flow;  (f)  (DCD(4C)l>)  flow,  (q)  (P|(ir|Ol  flow;  (»,)  l(7C)0)  flow,  (II  <  HC )  flow. 

Exp«r imental  condition*.  (*>:  R  "  2.54,  S  ►  41.2,  h  -  1.90  mm,  v  »  0.806  St . i  (b),  <c),  (d),  If),  (Mi  R  * 
0.213,  S  -  44.9,  h  -  1*82  mm,  v  •  0.789  St . »  (e),  (9):  R  -  0.282,  S  -  39.9,  h  «  1.93  mm,  v  -  0.784  St . t  (1), 
R  -  0.251.  S  ®  42.5.  h  -  1.87  mm.  v  -  0.792  St. 


At  this  stage  it  seems  worthwhile  to  attempt  a  summary  of  the  various 
flow  possibilities  correlated  to  the  flux,  the  quantity  that  was 
systematically  varied  throughout  the  experiment.  This  is  done  in  figure  13 
where  the  range  of  stability  of  the  various  flows  is  shown  as  a  function  of 
R.  To  keep  the  diagram  relatively  uncompl icated ,  flows  have  been  lumped 
together  in  classes,  which  correspond  to  the  sub  headings  of  §4.2,  and  an 
estimate  of  the  number  of  different  possible  flows  within  a  given  class  is 
indicated  on  the  right-hand  side  of  the  figure.  Of  course  figure  13  is  a 
gross  simplification  of  a  figure  that  should  be  drawn  in  a  five-dimensional 
parameter  space,  but  it  nevertheless  gives  a  rough  idea  of  the  various 
possibilities.  The  crinkled  lines  in  this  figure  indicate  what  I  shall  refer 
to  as  unsteady,  chaotic  motions. 


Eight-site  flows 

Flows  emanating  from  eight  sites  on  the  underside  of  the  plate  have 
already  been  encountered  in  figures  7  (i)  and  12  (d),  and,  as  indicated  by 
figure  13,  flows  in  this  class  were  observed  over  a  rather  wide  range  of  flow 
rates.  It  is  shown  in  figure  14  that  flows  with  (8-3)  drip  sites  and  (j)C 
streams,  for  3  =  0,...,8,  were  stable  possibilities  and  it  was  my  conviction 
while  carrying  out  the  experiment,  that  all  256  combinations  of  such  flows 
could  easily  be  realized.  I  did  not  try  to  delineate  the  ranges  of  stability 
for  each  of  the  various  kinds  of  eight-site  flows  though  there  were  evidently 
considerable  differences  between  flows.  For  example,  at  R  =  0.209,  S  =  47.4 
h  =  1.77  mm  and  v  =  0.799  St.,  all  combinations  of  flows  involving  zero  to 
six  drip  sites  seemed  possible,  but  I  was  unable  to  establish  stable  flows 
with  either  seven  or  eight  drip  sites.  On  the  other  hand,  when  the  flow  rate 
was  increased  slightly,  such  that  R  =  0.254,  S  =  41.2,  h  =  1.90  mm  and 
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variations  of  up  to  2%  in  the  flow  rate,  referred  to  in  §2: 


it  is  not  clear 


whether  the  unsteady  motions  arose  because  of  inherent  variations  in  the  flux 
or  whether  there  was,  in  fact,  an  indefinitely  persistent,  unsteady  motion 
corresponding  to  steady  operating  conditions  and  which,  in  the  experiments, 
eventually  lost  stability  to  a  more  stable  flow  merely  because  the  operating 
conditions  were  not  held  steady  enough. 

When  the  flow  rate  was  reduced  below  a  value  equivalent  to  R  ~  0.79  the 
streams  shown  in  figure  10  started  to  perform  a  complicated  'dance'  (see  the 
time-sequence  of  photographs  in  figure  12  (a)  -  (d),  and  a  sequence  under 
different  operating  in  (e)  -  ( i ) )  as  though  to  provide  room  for  an  eighth 
stream  to  form.  This  can  be  seen  clearly  in  figure  12  (h).  However  (cf. 
figure  12  ( i ) )  an  eighth  stream  did  not  develop  and  the  flow  settled  back 
towards  the  uniformly-spaced,  steady  (7C)  flow;  but  this  flow  was  not  stable 
...  .  A  typical  period  for  this  chaotic  type  of  motion  was  about  8  to  10  s, 
but  on  some  cycles  it  was  as  large  as  35  or  40  s.  The  unsteady  flow  shown  in 
figures  12  (a)  -  (c)  was  sustained  for  approximately  30  minutes  before 
collapsing  to  the  steady  (8C)  flow  shown  in  figure  12  (d).  The  unsteady  (7C) 
flow  shown  in  (e)  -  (i)  was  sustained  for  approximately  50  'cycles'  of  motion 
before  collapsing  to  a  steady  (8C)  flow. 

Thus,  as  indicated  above,  the  interpretation  of  these  observations  is 
somewhat  complicated,  but  it  is  my  view  that  an  unsteady,  chaotic  kind  of 
motion  would  have  arisen  and  persisted  indefinitely  if  the  operating 
conditions  had  been  held  more  nearly  steady.  We  shall  see  presently  that 
similar  kinds  of  unsteady  flows  arose  in  less  ambiguous  circumstances  in 
connection  with  other  flow  structures.  The  present  observations  indicate,  on 
the  other  hand,  the  importance  of  taking  into  account  'imperfections'  in 
maintaining  prescribed  operating  conditions  in  a  situation  in  which  there  are 
a  number  of  possible  flow  states. 
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FIGURE  12  ( ctd ) .  (e)-(i):  R  =  0.764 


FIGURE  19.  Some  unsteady  (T)  flows.  (a)  The  unsteady  (T)  flow  lost  stability  to  tills  steady  (7C)  flow  at 
S  =  13.2,  h  =  3.35  mm,  v  -  0.791  St.i  (b),  (c)  R  ■=  1.51,  S  =  12.9,  h  »  3.40  mm,  v  =  0.792  St.i 
S  =  12.7,  h  -  3.42  non,  u  =  0.791  St. 


to  R  ~  1-51#  the  unsteadiness  was  manifested  in  the  form  shown  in  figures  19 
(b),  (c)  where  a  'blob'  of  liquid  can  be  seen  to  have  developed  near  the  upper 
end  of  the  left-hand  side  of  the  sheet,  as  though  a  separate  stream  would  like 
to  form  there.  However,  this  'blob'  would  run  down  the  arm,  as  shown  in 
figure  19(c),  into  the  reservoir.  Then  the  left-hand  arm  might  pass  through  a 
phase  in  which  it  oscillated  back  and  forth  and  eventually  another  'blob' 
would  form.  This  chaotic  type  of  motion  seemed  to  be  quite  persistent,  being 
observed  for  over  fifteen  minutes  before  the  operating  conditions  were 
changed.  A  further  reduction  of  the  flow  rate  resulted  in  these  unsteady  (T) 
flows  losing  stability  to  the  steady  (7C)  flow  shown  in  figure  19(a). 

(C  +  T)  flows 

A  class  of  flows  consisting  of  one  continuous  get  of  liquid  and  a  'T' 
section,  similar  in  form  to  that  described  under  the  class  of  (T)  flows,  was 
observed  at  flow  rates  equivalent  to  R  aproximately  in  the  range  [1.50, 
2.16].  Some  examples  from  the  branches  of  steady  (C+T)  flows  are  given  in 
figure  20:  (b)  shows  a  (CT)  flow  and  (c)  is  a  (TC)  flow  at  the  same  operating 

conditions  as  for  (b);  (d)  and  (e)  show  (CT)  flows  at  larger  flow  rates;  and 

this  branch  of  (CT)  flows  lost  stability  at  R  =  2.16,  S  =  10.1,  h  =  3.84  mm, 
v  =  0.795  St.  to  the  steady  (T)  flow  shown  in  (f).  The  flow  depicted  in 
figure  20(a)  at  R  =  1.57,  S  =  12.4,  h  =  3.46  mm,  v  =  0.797  st.  was  near  the 
limit  of  stability  for  steady  (CT)  flows,  and  a  further  small  reduction  in  the 
flow  rate  led  to  the  development  of  an  unsteady  (CT)  flow,  examples  of  which 
are  given  in  figure  21.  The  unsteadiness  in  (d)  and  (e)  took  the  form  of  the 
left-hand  arm  of  the  'T'  section  oscillating  back  and  forth  in  the  plane  of 
the  'T',  with  (d)  showing  the  flow  in  a  relatively  quiescent  phase  of  its 
motion  and  (e)  in  an  'active'  phase  (as  evidenced  by  the  reflection  of  the 
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left-hand  arm  of  the  'T'  section  in  the  free  surface  of  the  reservoir).  Still 
further  reductions  of  the  flow  rate  led  to  a  more  vigorous  form  of  the 
unsteadiness  as  indicated  in  (b)  and  (c),  both  of  which  were  taken  under  the 
same  operating  conditions,  and  in  (a)  where  a  liquid  'blob',  similar  to  that 
observed  in  the  unsteady  (T)  flows,  can  be  seen  on  the  left-hand  arm  of  the 
'T'  section.  (Note  that,  in  (a),  (b)  and  (c)  the  two  arms  at  the  edges  of  the 
'T*  section  came  together  near  the  level  of  the  free  surface  in  the 
reservoir.)  The  branch  of  unsteady  (CT)  flows  lost  stability  to  a  steady  (7C) 
flow  at  R  =  1.50,  S  =  12.8,  h  =  3.41  mm,  v  =  0.797  St.,  corresponding  to  a 
flow  rate  just  below  that  which  obtained  for  figure  21(a). 

( 2C  +  T)  flows 

Examples  of  steady  flows  from  this  class  are  shown  in  figures  22(b)  - 
(d).  The  photograph  of  figure  22(a)  is  an  oblique  view,  almost  a  side 
elevation,  showing  how  the  two  arms  of  the  'T*  section  came  together  to  form  a 
kind  of  fan  shape  in  the  plane  orthogonal  to  the  plane  of  the  triangular  'T' 
section. 

Further  examples  from  the  branch  of  steady  (CTC)  flows  are  given  in 
figure  23.  This  branch  lost  stability  with  increasing  flow  rate  at  R  =  2.06, 
S  =  10.4,  h  =  3.79  mm,  v  =  0.798  St.  to  the  steady  (CT)  flow  shown  in  figure 
23(c).  The  branch  lost  stability  with  decreasing  flow  rate  at  R  =  1.64,  S  = 
12.1,  h  =  3.51  mm,  v  =  0.798  St.  to  a  branch  of  unsteady  (CTC)  flows,  an 
example  of  which  is  given  in  figure  24.  Shown  there  is  a  sequence  of  photo¬ 
graphs,  under  fixed  operating  conditions  for  the  flow,  which  illustrate  the 
rather  violent  nature  of  the  gyrations  of  the  central  'T'  section  of  the 
pattern,  while  the  'C'  streams  at  either  side  remained  essentially  unperturbed 
by  the  proceedings.  This  branch  of  unsteady  (CTC)  flows  collapsed  at  R  = 
1.55,  S  =  12.5,  h  =  3.45  mm,  v  =  0.798  St.  to  an  unsteady  (CCTC)  pattern,  and 


FIGURE  23.  Examples  fro*  the  branch  of  steady  (CTC)  flows. 
1.94, 
flow 


(a)  R  -  1.71, 


11.7,  h 


*  -  -  - r  w  -  . .  -  _>,  v  -  0.799  St.l  (b)  R  • 

S  -  TO. 7,  h  -  3.72  mm,  w  -  0. 7o9  St.,  (c)  the  steady  (CTC)  flows  lost  stability  to  this  steady  (CT) 
at  R  •  2,06,  fl  -  10.4,  h  •  3.79  ■»,  v  -  0.790  Bt . 
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this  changed  spontaneously  after  about  2  V2  min  to  a  steady  ( 7 C )  flow. 

Flows  of  the  (CCT)  kind  had  similar  properties  to  those  just  described. 
The  branch  of  steady  (CCT)  flows,  typified  by  the  photographs  of  figure  25, 
lost  stability  with  increasing  flow  rate  at  R  =  1.91,  S  =  10.9,  h  =  3.69  mm, 
v  =  0.796  St.  to  the  steady  (T)  flow  shown  in  figure  25(c).  As  the  flow  rate 
was  decreased  the  steady  (CCT)  flow  developed  into  an  unsteady  motion  at  R  = 
1.59,  S  =  12.3,  h  =  3.47  mm,  v  =  0.795  St.  (see  figure  (26(e)).  Another 
example  of  an  unsteady  (CCT)  flow  is  shown  in  figures  26(b)  -  (d)  in  a 
sequence  of  photographs  taken  at  a  given  operating  state.  The  branch  of 
unsteady  (CCT)  flows  lost  stability  at  R  =  1.50,  S  =  12.9,  h  =  3.40  mm, 
v  =  0.794  St.  to  the  steady  (7c)  flow  shown  in  figure  26(a). 

(3C  +  T)  flows 

Flows  involving  three  continuous  jets  of  liquid  and  one  ’T’  section  were 
found  to  be  stable  at  flow  rates  equivalent  to  R  lying  between  about  1.5  and 
1.8.  Examples  of  steady  (CCTC)  flows  are  given  in  figures  27(a),  (b),  (d), 
and  an  example  of  a  steady  (CTCC)  flow  is  shown  in  figure  27(c).  The  branch 
of  steady  (CCTC)  flows  lost  stability  with  increasing  flow  rate  at  R  =  1.86, 

S  =  11.0,  h  =  3.68  mm,  v  =  0.802  St.  to  the  steady  (TC)  flow  shown  in  figure 
27(e).  On  decreasing  the  flow  rate,  the  branch  developed  into  an  unsteady, 
chaotic  type  of  motion  at  parameter  values  near  R  =  1.60,  S  =  12.2,  h  =  3.49 
mm,  v  =  0.801  St.,  some  examples  of  the  flow  at  these  conditions  being  given 
in  figures  28(e),  (f).  Most  of  the  unsteadiness  in  this  flow  was  associated 
with  pulsations  of  the  triangular  section  or  with  sideways  movem"nts  of  the 
'leg'  beneath  the  triangular  section,  both  features  of  which  are  evident  in 
the  temporal  sequence  of  photographs  shown  in  figures  28(b)  -  (d).  In  spite 
of  the  vigorous  motions  exhibited  by  this  flow,  the  pattern,  nevertheless, 
preserved  its  overall  coherence,  returning  from  time  to  time  to  a  nearly 


dy  ( CCT )  flows.  (a)  Stability  was  lost  to  this  ( 7C )  flow  at  P.  *•  1.59,  S  =  12.3,  h  *  3 
St.;  (b),  (c),  (d)  are  a  sequence  In  time  of  an  unsteady  (CCT)  flow  at  R  -  1.54,  S  =  1 
=  0.792  St.;  (e)  unsteady  (CCT)  flow  at  R  *  1.50,  S  =  12.9,  h  =  3.40  inm,  «  -  0.794  St. 


Examples  from  the  branch  of  ((V^TJQ)  flows  are  given  in  figure  37. 

This  branch  lost  stability,  with  decreasing  flow  rate,  to  the  steady  (T)  flow 
shown  in  figure  37(a)  and,  with  increasing  flow  rate  at  R  =  4.50,  S  =  6.16, 
h  =  3.58  mm,  v  =  0.796  St.,  to  a  ( QQ )  flow.  The  sensitivity  of  these  flows  to 
parameters  other  than  R  is  indicated  by  a  comparison  between  figures  37(c) 
and  (d),  taken  at  nearly  the  same  values  of  R  but  with  quite  different  fluid 
viscosities . 


Lower-faa-  flows 

As  the  flow  rate  was  increased  a  stage  was  i .  a  is  s-s't  i  bed  above, 

when  the  liquid  curtain  completely  spanned  *'■  :  >:  u:.f  inker,  sheet, 

disconnecting  the  region  in  front  of  the  .r  •  a  .  •  :  :  •  •  When  these 

circumstances  prevailed,  the  cavity  S,eh.  •  :  r  •  •  •••...•  ,  »  •  wa.  pounded  by 

the  walls  of  the  reservoir,  the  plate,  ; :  *  *  .  t  ::«•  reservoir  and 

the  curtain  itself,  was  vented  to  the  a '  »«*  si  •  •  ■  t  ••  t  ..re  *  :  e  j  t  assure  on 

either  side  of  the  curtain.  This  was  done  !  y  !■*  i.fu.g  a  piece  of  1  cm-bore 
cupper  tubing  in  the  form  of  a  s>  r  fiat  it  hugged  the  front  and  rear  walls 

and  the  bottom  of  the  reset  voi  i  .  Htgr.  :  well  away  from  the  zone  at  which  the 
liquid  sheet  penetrate!  t  he  free  surface  of  the  reservoir,  this  tube  had 
negligible  effect  on  t  he  tl  -w  fie’..',  ir.  the  reservoir. 

The  view  looking  t  war  is  f  *'e  er.i  of  the  plate  revealed  very  little  for 
this  class  if  flows,  ar  i  «  .’all  tiy  to  ■ -Pa  rue  ter  i  ze  them  from  the  i  r 
appjearance  from  t  fie  sid»-,  isir.-i  *  ;  'Sil  io*.  <  ;  t  tie  a  1 1  achmen  f  line  of  the 

free  surface  on  the  undersoil-  of  h'i-  j  late  as  a  dist  inguishing  feature  of  the 
flow.  Thus,  the  term  *  lowe i - f a-  e  f  Paw1  is  used  to  t  i-for  to  those  flows  for 
which  a  falling,  unbroken  sheet  of  liquid  completely  spanned  the  reservoir  and 


for  which  the  line  of  attachment  of  the  free  sut  I  <i'-f  with  t  tie  plate  was  on  the 


r 
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Quadrant  flows 

The  term  'quadrant  flow'  has  been  used  to  indicate  those  flows  in  which  a 
free  boundary  of  the  liquid  curtain  attached  to  a  side  wall  of  the  reservoir, 
leaving  a  'hole'  in  the  curtain  shaped  somewhat  like  a  quadrant  of  a  circle. 
Examples  of  such  flows  have  already  been  given  under  the  heading  of 
'partially-attached  flows’  and  three  further  examples  showing  a  (QfV^T)),  a 
(QQ)  and  a  ((V^TJQ)  flow,  all  established  at  the  same  operating  conditions, 
are  given  in  figure  35.  All  the  flows  observed  within  this  class  were  steady. 

The  structure  of  the  branch  of  (QQ)  flows  is  illustrated  in  the 
collection  of  photographs  in  figure  36.  As  the  flow  rate  was  increased,  the 
attachment  points  of  the  'quadrants'  moved  down  the  side  walls  leaving 
increasingly  smaller  'holes'  in  the  liquid  curtain  -  only  very  small  holes 
remained  at  the  conditions  for  figure  36(g)  -  until  eventually  a  complete, 
unbroken  sheet  had  formed  at  a  flow  rate  corresponding  approximately  to  R  = 
6.33,  S  =  5.27,  h  =  5.31  mm,  v  =  0.755  St.  On  decreasing  the  flow  rate  the 
attachment  points  moved  up  the  side  walls  of  the  reservoir  until,  with  R  = 
4.18,  S  =  6.93,  h  =  4.63  mm,  v  =  0.756  St.,  the  attachment  point  of  the  right- 
hand  quadrant  had  reached  the  plate  (cf.  figure  36(a)).  It  is  interesting 
that  a  reduction  in  flow  rate  of  about  10%  accounted  for  the  relatively  small 
changes  evident  between  (c)  and  (b)  of  figure  36,  whereas  the  flow  rates  for 
(b)  and  (a)  differed  by  only  about  1.5%.  Some  side  views  of  (QQ)  flows  are 
shown  in  figures  36(d)  and  (e):  in  (d)  the  'quadrant'  attached  to  the  near 
wall  at  about  half  way  between  the  plate  and  the  level  in  the  reservoir;  in 
(e)  the  attachment  to  the  near  wall  was  very  close  to  the  junction  with  the 
free  surface. 
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flow  that  was  observed  (see  figure  33)  was  an  unsteady  flow  and,  in  spite  of 
apparently  having  a  small  'domain  of  attraction',  seemed  to  be  quite  a  stable 
motion  as  it  persisted  for  over  15  minutes,  at  which  time  a  different  flow  was 
induced . 


Partially-attached  flows 

As  the  flow  was  increased  in  this  experiment  there  was  a  tendency  for  the 
fluid  more  and  more  to  fall  as  a  curtain,  or  sheet,  of  liquid.  This  was 
apparent  with  the  class  of  (T)  flows  described  above  and  is  further 
exemplified  in  the  classes  of  motions  to  be  described.  Indeed,  once  the  flow 
rate  was  sufficiently  large  it  was  possible  to  establish  a  flow  with  the 
'sheet'  attached  to  one  wall  of  the  reservoir,  but  not  necessarily  to  both 
walls,  and  at  even  larger  flow  rates  only  those  flows  with  the  'sheet' 
attached  to  both  walls  were  possible.  States  with  the  'sheet'  attached  to  one 
wall  only  will  be  referred  to  as  partially-attached  flows,  some  examples  of 
which  are  given  in  figure  34.  The  liquid  'curtain'  in  figure  34(a)  was 
attached  to  the  right-hand  wall  and  a  small  'quadrant'  of  the  sheet  was  left 
open  at  the  corner  near  the  free  surface  in  the  reservoir  at  the  left-hand 
wall;  flows  of  this  kind  will  be  referred  to  as  (QA)  flows.  The  photograph  in 
figure  34(b),  taken  at  the  same  parameter  settings  as  for  (a),  shows  the 
liquid  curtain  attached  to  the  left-hand  wall,  with  the  right-hand  boundary  of 
the  curtain  having  the  appearance  of  the  right-hand  half  of  a  (T)  flow.  Such 
flows  will  therefore  be  denoted  as  (AtV^T))  flows. 

Some  side  views  of  these  kinds  of  flows  are  shown  in  figures  34(c)  -  (e), 
from  which  it  is  apparent  that  the  attachment  line  with  the  plate  was  on  the 
underside,  well  behind  the  lower  corner  of  the  plate.  Partially-attached 
flows  were  found  over  a  rather  wide  range  of  flow  rates,  as  indicated  in 
figure  13.  They  were  all  steady  motions. 
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in  an  unsteady  form  and  there  appeared  not  to  be  an  associated  branch  of 


steady  flows  of  this  pattern.  Examples  of  the  ( TCT)  flows  are  given  in  figure 
31.  At  each  operating  condition  a  pair  of  photographs  is  included  to  indicate 
the  nature  of  the  unsteadiness.  Thus,  not  only  did  the  triangular  section 
pulsate  as  in  previous  flows,  but  the  central  jet  also  moved  to  and  fro  in 
plane  of  the  falling  liquid.  On  increasing  the  flow  rate  this  branch  of  flows 
lost  stability  at  R  =  1.90,  S  =  11.0,  h  =  3.67  mm,  v  =  0.793  St.  to  the 
steady  (T)  flow  shown  in  figure  31(1.).  On  reducing  the  flow  rate  the  (TCT) 
flows  lost  stability  at  R  =  1.54,  S  =  12.7,  h  =  3.43  mm,  v  =  0.792  St.  to  the 
(CCT)  flow  shown  in  figure  31(a),  which  was  also  an  unsteady  motion. 

(2T  ■+  2C)  flows 

As  in  the  case  of  (TCT)  flows,  the  present  class  of  flows  were  observed 
only  as  unsteady  motions.  One  such  flow,  a  (TCTC)  flow,  is  illustrated  in 
figure  32.  The  unsteadiness  in  this  flow,  which  is  apparent  by  comparisons 
between  the  pairs  of  photographs  at  the  same  operating  conditions,  took  the 
form  of  a  pulsation  of  the  two  'T'  sections  on  a  time  scale  of  about  1  s  or 
so,  and  of  a  wandering  to-and-fro  of  the  1 C*  stream  between  the  two  ’T' 
sections  on  a  time  scale  of  around  10  to  15  s;  the  right-hand  jet  remained 
relatively  unperturbed.  This  branch  of  flows  lost  stability,  with  increasing 
flow  rate,  at  R  =  1.88,  S  =  11.1,  h  =  3.67  mm,  v  =  0.795  St.  to  the  steady 
(TC)  flow  shown  in  figure  32(g).  Unfortunately  I  did  not  determine  how  the 
branch  of  (TCTC)  flows  lost  stability  with  decreasing  flow  rate,  but  it 
appeared  as  though  the  branch  of  (CTCT)  flows  lost  stability  to  a  steady  (CCT) 
pattern  at  values  of  the  parameters  of  approximately  R  =  1.72,  S  =  11.8,  h  = 
3.55  mm,  v  =  0.792  St. 

On  one  occasion  I  was  able  to  establish  a  (TCCT)  flow  but  subsequently  1 


was  unable  to  manipulate  the  conditions  to  reestablish  this  flow.  The  (TCCT) 


steady  (CCTC)  flow.  This  branch  of  unsteady  (CCTC)  flows  lost  stability  at 
K  -  1.48,  S  =  12.9,  h  =  3.40  mm,  v  =  0.799  St.  to  the  steady  (7C)  flow  shown 
in  f igure  28(a). 

All  attempts  to  produce  a  pattern  of  the  form  (TCCC)  or  (CCCT)  were 
unsuccessful . 


( 4C  +  T)  flows 

The  only  flow  patterns  within  this  class  of  flows  that  I  was  able  to 
establish  were  those  of  the  form  (CCTCC),  some  examples  of  the  steady  version 
of  this  flow  being  given  in  figure  29.  This  branch  of  flows  lost  stability, 
with  increasing  flow  rate,  at  R  =  1.91,  S  =  11.1,  h  =  3.66  mm,  v  =  0.787St. 
to  the  steady  (T)  flow  shown  in  figure  29(d).  On  decreasing  the  flow  rate  the 
flow  developed  an  unsteady,  chaotic  type  of  motion  at  values  of  the  parameters 
near  R  =  1.62,  S  =  12.4,  h  =  3.47  mm,  v  =  0.787  St.  The  nature  of  the 
unsteadiness  of  the  pattern  at  these  conditions  is  indicated  in  figures  30(f), 
(g),  and  a  more  vigorous  form  of  the  unsteadiness  is  illustrated  by  the  photo¬ 
graphs  of  fiyures  30(c),  (d)  and  (e).  This  sequence  shows  how  the  'leg'  of 
the  triangular  section  might  undergo  sideways  oscillations,  or  how  the 
triangular  section  itself  might  pulse  'up  and  down'.  This  branch  of  unsteady 
(CCTCC)  flows  lost  stability  at  R  =  1.53,  S  =  12.6,  h  =  3.43  mm,  V  =  0.798 
St.  to  the  steady  (7C)  flow  shown  in  figure  30(b).  An  early  stage  in  the 
transition  to  the  steady  (7C)  flow  is  shown  in  figure  30(a). 

( TCT )  flow 

At  flow  rates  equivalent  to  R  lyiny  between  roughly  1.6  and  1.9  there 
was  found  a  branch  of  unsteady,  chaotic  type  of  flows  consisting  of  two  'T' 


portions  separated  by  a  continuous  jet  of  liquid.  These  flows  were  seen  only 


lower  face  of  the  plate.  One  such  branch  of  flows,  to  be  denoted  (LF1)  flows, 


is  illustrated  in  figure  38.  This  branch  of  flows  arose  naturally  from  the 
class  of  'quadrant  flows'  at  around  a  flow  rate  equivalent  to  R  =  6.3.  Once 
established,  however,  the  branch  of  ( LF 1 )  flows  was  stable  down  to  a  flow  rate 
approximately  equivalent  to  R  =  2.59,  S  =  9.10,  h  =  4.04  mm,  v  =  0.784  St. 
Unfortunately  I  do  not  have  any  photographs  of  side  views  of  ( LF 1 )  flows  at 
values  of  R  smaller  than  that  applying  to  figure  38(a),  namely  R  =  3.50, 
but  the  appearance  of  these  flows  was  not  very  different  from  the  photographs 
shown  in  figures  34(c)  and  (d)  for  'partially-attached  flows'. 

Being  focussed  on  a  region  close  to  the  near  wall  of  the  reservoir, 
photographs  such  as  those  shown  in  figure  38  give  on  incomplete  picture  of  the 
flow.  For  these  particular  flows  the  liquid  curtain  was  actually  bowed,  and 
near  the  side  walls  the  attachment  point  of  the  meniscus  was  much  closer  to 
the  lower  corner  of  the  plate  then  it  was  in  the  centre  of  the  channel.  This 
is  indicated  schematically  in  figure  39(a).  In  figures  38(a),  (b)  the 
'bowing'  of  the  curtain  was  not  very  large,  but  for  the  conditions  obtaining 
for  figures  38(c)  -  (f)  it  was  appreciable.  For  example,  in  figure  38(d)  the 
attachment  oint  at  the  centre  of  the  channel  was  displaced  about  0.6  to  0.7a 
further  from  the  corner  than  was  the  attachment  point  near  the  well.  It  is 
this  effect  that  accounts  for  the  out-of-focus  structure  to  the  rear  of  the 
sharp,  focussed  lines  of  the  meniscus  on  the  near  wall  of  the  reservoir. 

Also,  the  dark  patch  below  the  free  surface  in  the  reservoir  corresponds  to 
the  indentation  of  the  free  surface  near  the  entry  region  for  the  curtain,  and 
the  lateral  extent  of  this  patch  (note  that,  being  out  of  focus,  it  too 
appears  blurred)  is  indicative  of  the  amount  of  bowing  in  the  curtain. 

The  attachment  point  of  the  free  surface  in  the  zone  near  the  walls  of 
the  reservoir  was  approximately  0.57a  from  the  corner  of  the  plate  when 
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FlCUhfc  _S  i.  Schematic  lej-rtientanon  of  the  attachment  line  on  the  underside  of  the  plate  for  (LFl  )  and  (LF2)  flows. 

(a)  Typical  jrofiie  for  (LFl )  flows;  (b)  showing  the  transition  from  an  ( LF2 )  profile  to  an  (LFl)  profile 
(fc),  (r_  j  ,  ;  i)  s;,ow  yi.en.at  leal  ly  how  the  (LF2)  profiles  changed  with  increasing  flow  rate. 
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R  =  3.50,  as  in  figure  38(a).  With  R  =  4.67,  as  for  figure  38(b),  this 
attachment  point  had  moved  to  a  distance  of  approximately  0.75a  from  the 
corner  and,  with  R  =  6.61,  the  distance  had  increased  to  about  0.98a  (cf. 
figure  38(c)).  Near  ttie  stability  boundary  for  the  branch  of  (LF1)  flows,  the 
attachment  point  close  to  the  walls  of  the  reservoir  (see  figure  38(f))  was 
about  1.15a  from  the  corner  of  the  plate.  A  small  increase  in  the  flow  rate 
from  the  conditions  obtaining  for  figure  38(f)  brought  about  a  major  change  in 
the  flow  structure,  as  indicated  by  figure  38(g),  with  the  attachment  point 
moving  forward  to  a  position  only  about  0.27a  from  the  lower  corner  of  the 
plate.  This  new  branch  of  flows  will  be  referred  to  as  the  (LF2)  branch. 

Examples  of  the  ( LF2 )  branch  of  flows  are  given  in  figure  40,  the  photo¬ 
graphs  being  ordered  by  increasing  flow  rate  (or  equivalently  by  increasing 
R,  since  v  was  nearly  constant  for  the  sequence).  The  (LF1)  branch  of 
flows  lost  stability  to  the  ( LF2 )  branch  at  R  =  10.1,  S  =  3.93,  h  =  6.15  mm, 
v  =  0.745  St.,  corresponding  to  operating  conditions  lying  between  (c)  and  (d) 
in  the  ordering  of  figure  40.  As  indicated  above,  the  tree  surface  attached 
much  closer  to  the  corner  of  the  plate  for  the  (LF2)  flow  shown  in  figure 
38(g)  than  was  the  case  for  the  (LF1)  flow  of  figure  38(f).  Another 
difference  between  ( LF 1 )  and  (LF2)  flows  was  the  shape  of  the  attachment  line 
on  the  underside  of  the  plate,  a  difference  which  is  represented  schematically 
in  figure  39.  Thus,  the  curvature  of  the  profile  of  the  attachment  line  for 
(LF1)  flows  was  one-signed,  whereas  the  profile  for  ( LF2 )  flows  had  five  local 
extrema  (on  an  appropriate  open  interval).  This  difference  accounts  for  an 
apparent  slight  'fuzziness'  of  the  free  surface  attaching  to  the  underside  of 
the  plate  in  figure  38(g),  and  in  the  photographs  of  figure  40,  compared  with 
the  sharp  outline  evident  in  figure  38(f),  sav .  ( Aside .  My  notes  and  my 


recollection  of  the  experiment  are  consistent  with  the  above  description,  but 
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the  texture  of  the  curtain  in  figure  38(f)  raises  the  possibility  that  the 
attachment  profile  of  the  (LF1)  flow  may  have  developed  more  'structure'  than 
I  have  indicated  above,  as  the  stability  boundary  was  approached,  'anticipat¬ 
ing'  the  transition  to  an  ( LF2 )  flow.) 

After  the  (LF2)  flow  shown  in  figure  38(g)  had  been  established, 
reductions  in  the  flow  rate  brought  about  changes  typified  by  figures  40(c) 
and  (b) ,  and  represented  schematically  in  figure  39  with  regard  to  the 
attachment  profile.  When  the  flow  rate  was  adjusted  slightly  below  that 
obtaining  for  figure  40(b)  the  liquid  curtain  first  developed  an  extra  pair  of 
corrugations,  as  indicated  by  the  dashed  'attachment  profile',  and  then  the 
motion  evolved  into  the  (LF1 )  flow  shown  in  figure  40(a). 

Changes  in  the  structure  of  the  branch  of  (LF2)  flows  with  increasing 
flow  rate  are  exemplified  in  figure  40  (and  cf.  figure  39(b)  -  (d)).  In 
particular,  the  attachment  line  on  the  underside  of  the  plate  moved  closer  and 
closer  to  the  (lower)  corner  of  the  plate  and  the  liquid  curtain  became  more 
nearly  vertical,  until  at  the  conditions  obtaining  for  figure  40 (k)  the 
attachment  line  was  located  very  nearly  at  the  corner  of  the  plate.  Then,  a 
slight  increase,  of  less  than  2%,  in  the  flow  rate  resulted  in  a  dramatic 
change  in  flow  structure  to  that  shown  in  figure  40(1),  with  the  meniscus  now 
attached  to  the  (nearly)  vertical  face  of  the  plate.  This  new  class  of  flows 
will  be  referred  to  as  (VF)  flows. 

Vertical-face  flows 

The  transition  from  the  branch  of  (LF2)  flows  to  the  (VF)  flow  shown  in 
figure  40(1)  is  illustrated  in  figure  41.  Here  the  end  view  of  the  flow  shows 
the  profile  of  the  meniscus  line  on  the  (nearly)  vertical  face  of  the  plate. 
When  the  branch  of  (LF2)  flows  lost  stability,  the  meniscus  line  on  the  plate 
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moved  onto  the  end  surface  of  the  plate  and  developed  a  number  of 
corrugations,  as  shown  in  figure  41(a).  But  this  corrugated  structure 
evidently  was  not  an  equilibrium  configuration,  and  the  corrugation  near  the 
right-hand  wall  of  the  reservoir  moved  slowly  towards  this  wall.  On  reaching 
the  sharp  gradient  in  the  meniscus  line  near  the  wall  of  the  reservoir,  the 
corrugation  disappeared,  leaving  a  pattern  of  the  form  shown  in  figure  40(b). 
Surface-tension  forces  precluded  the  flow  of  figure  40(b)  from  being  an 
equilibrium  state,  and  slowly  the  corrugation  near  the  left-hand  wall  moved 
sideways  and  disappeared  in  the  same  manner  as  above,  leaving  the  pattern 
shown  in  figure  41(d).  The  remaining  two  corrugations  then  disappeard,  not  by 
moving  to  the  edges  of  the  reservoir  but  by  diminishing  in  amplitude  while 
remaining  fixed  horizontally,  to  leave  the  equilibrium  pattern  shown  in  figure 
41(e).  The  latter  flow  will  be  referred  to  as  a  (VF(0))  flow,  the  zero 
indicating  that  the  attachment  line  on  the  end  of  the  plate  was  not 
corrugated.  (Note :  the  vertical  tube,  somewhat  out  of  focus,  in  the  fore¬ 
ground  on  the  right-hand  side  of  the  photographs  is  part  of  the  'U-tube'  used 
to  equalize  the  pressure  on  the  two  sides  of  the  liquid  curtain.) 

Further  examples  of  the  VF(0)  branch  of  flows  are  given  in  figure  42. 

The  photographs  shown  in  figures  42(b),  (c)  and  (d)  are  different  views  of  the 
same  flow:  (b)  is  a  side  view  showing  how  the  meniscus  attached  to  the  plate 
and  (d)  indicates  the  shape  of  the  meniscus  line  on  the  end  of  the  plate;  the 
way  the  meniscus  attached  to  the  plate  in  the  zone  near  the  side  walls  of  the 
reservoir  is  indicated  in  (c),  where  it  is  seen  that  the  attachment  line  of 
the  free  surface  passed  from  the  end  face  of  the  plate  to  its  lower  face.  As 
the  flow  rate  was  increased  tls  transition  point  moved  closer  to  the  side 
walls  (cf.  figures  42(e),  (g),  (  i )  ) . 
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Whereas  for  (LF2)  flows  the  falling  curtain  was  nearly  vertical  (cf. 
figure  40),  it  extended,  in  the  case  of  (VF)  flows,  well  beyond  the  end  of  the 
plate.  The  location  of  the  meniscus  on  the  side  wall  of  the  reservoir  is 
evident  in  figure  42(b),  and  the  spatial  location  of  the  curtain  near  the 
central  region  of  the  channel  is  indicated  by  the  dark  band  at  the  front  of 
the  photograph  (cf.  also  figure  42(f)  taken  under  different  lighting 
conditions).  Thus,  the  curtain  was  quite  strongly  curved  in  the  regions  near 
the  walls  of  the  reservoir. 

At  the  smaller  flow  rates  associated  with  this  branch  of  flows  the 
meniscus  line  on  the  end  of  the  plate  rose  to  a  maximum  height  and  then 
dropped  slightly  to  a  local  minimum  (cf.  figure  42(a),  where  there  is  also 
evidence,  on  the  left-hand  half  of  the  plate,  of  a  second  pair  of  local 
extrema).  As  the  flow  rate  was  increased  this  local  minimum  became  less 
pronounced  and  in  figure  42 ( i )  is  nearly  indiscernible.  The  flow  shown  in 
figure  42 ( i )  was  at  the  maximum  flow  rate  achievable  with  the  pump  used  for 
these  experiments. 

Following  the  loss  of  stability  of  the  (LF2)  flows  and  the  formation  of  a 
( VF ( 0 ) )  flow,  the  flow  rate  could  then  be  reduced  with  the  (VF(0))  flow  main¬ 
taining  stability.  In  one  experiment,  for  example,  the  (VF(0))  flow  formed 
at  R  =  28.0,  S  =  1.95,  h  =  8.74  mm,  v  =  0.758  St.,  and  this  branch  of  flows 
was  stable  on  reduction  of  the  flow  rate  until  R  =  27.0,  S  =  2.01,  h  =  8.59 
mm,  v  =  0.752  St,  at  which  stage  a  different  kind  of  (VF)  motion  developed. 
Another  realization  of  the  same  sequence  is  illustrated  in  figures  42  and  43: 
the  flow  shown  in  figure  42(a)  was  near  the  stability  margin  for  iVF(O))  flows 
and  a  small  reduction  in  flow  rate  led  to  the  unsteady  ( VF ( 0,1))  flow  shown  in 
figure  43.  The  latter  flow  is  so  designated  because,  as  shown  in  figure 
43(a),  at  a  certain  phase  of  the  motion  the  meniscus  line  on  the  end  plate  was 
free  of  corrugations,  whereas  at  a  later  phase  the  meniscus  line  developed  a 
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single  corrugation  near  the  right-hand  wall  of  the  ieservoir  (see  figure 
43(b)).  This  corrugation  then  moved  towards  the  side  wall,  as  shown  in  figure 
43(c),  eventually  losing  its  identity  and  allowing  the  pattern  shown  in  (a)  to 
reestablish  itself.  The  overall  motion  was  nearly  periodic  in  time,  on  a  time 
scale  of  about  8  s. 

As  the  flow  rate  was  reduced  the  (VF(0,1))  flow  persisted  until  the 
conditions  were  such  that  R  =  26.0,  S  =  2.06,  h  =  8.50  mm,  v  =  0.754  St.,  at 
which  stage  the  flow  illustrated  in  figure  44  developed.  This,  too,  was  a 
near ly-per iodic  motion  involving  either  five  or  six  corrugations  and  will  be 
referred  to  as  a  (VF(5,6))  flow.  That  phase  of  the  motion  where  there  were 
five  corrugations  is  illustrated  in  figure  44(a)  and  the  phase  involving  six 
corrugations  is  illustrated  in  figure  44(b).  It  is  seen  that  the  sixth 
corrugation  formed  at  the  right-hand  side  of  the  pattern.  Once  it  had  formed, 
the  new  corrugation  moved  slowly  towards  the  right-hand  wall  of  the  reservoir, 
losing  its  identity  when  it  reached  the  sharp  gradient  of  the  meniscus  line 
near  the  wall.  The  approximate  period  of  this  motion  was  about  40  s  and  it 
was  observed  for  approximately  16  minutes  before  the  flow  rate  was  changed. 
(Note  that  the  corrugating  of  the  free  surface  increases  the  mean  curvature, 
enabling  it  to  support  a  greater  load  than  would  otherwise  be  possible.)  The 
presence  of  the  corrugations  caused  the  uppe-  and  lower  surfaces  of  the  liquid 
curtain  to  become  wrinkled,  some  of  which  structure  can  be  seen  in  figure  44 
from  reflections  from  the  free  surface  (and  cf.  figure  42  where  the  free 
surface  was  much  less  wrinkled). 

A  small  reduction  in  the  flow  rate  from  that  obtaining  for  figure  44, 
giving  operating  conditions  of  R  =  25.6,  S  =  2.09,  h  =  8.43  mm,  v  =  0.751 
St.,  initiated  a  change  of  flow  structure  to  the  (VF(6,5))  flew  illustrated  in 
figure  45.  As  for  the  (VK(5,6))  flow  this  motion  also  was  nearly  periodic  in 
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time,  but  the  approximate  period  was  much  smaller,  taking  a  value  of  about 
4s.  In  addition  the  sixth  corrugation  now  formed  on  the  left-hand  side  of  the 
pattern  (cf.  figures  45(a)  and  (b))  and  proceded  leftwards  (cf.  figure  45(c)), 
losing  its  identity  near  the  left-hand  wall.  The  flow  shown  in  figure  45  was 
observed  for  20  minutes  before  the  flow  rate  was  reduced.  Further  small 
reductions  in  the  flow  rate  had  little  effect  on  the  general  features  of  the 
flow  illustrated  in  figure  45,  except  that  the  'period*  of  the  motion 
increased  slightly,  reaching  a  value  of  about  6  s  when  R  =  24.6,  S  «  2.18, 
h  =  8.25  Bin,  v  =  0.741  St.  Then  a  small  reduction  in  flow  rate,  to  the 
conditions  R=  24.3,  S  =  2.21,  h  *  8.21  mm,  v  *  0.741  St.,  increased  the 
'period'  of  the  motion  to  a  value  of  about  40  s.  With  a  further  small 
reduction  in  flow  rate,  such  that  R  =  23.9,  S  *  2.23,  h  =  8.16  mm,  v  =  0.740 
St.,  the  'period*  increased  to  around  100  s  and  the  sixth  corrugation  seemed 
to  have  considerable  difficulty  in  forming:  the  meniscus  line  'dimpled'  as 

though  the  corrugation  would  form,  but  it  did  not  and  the  line  straightened 
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out  again;  however,  after  about  30  to  40  such  attempts  the  sixth  corrugation 
would  finally  become  established  and  then  it  would  translate  slowly  towards 
the  left-hand  wall. 

When  the  flow  rate  was  then  reduced  to  a  value  such  that  R  -  23.7,  S  = 
2.25,  h  *  8.12  mm,  v  »  0.738  St.,  a  qualitative  change  occurred  in  the  flow 
structure,  as  indicated  in  figure  46.  In  this  flow  a  new  corrugation  would 
develop  near  one  wall  and  the  corrugation  near  the  opposite  wall  would  'slip 
away'.  So,  for  example,  a  new  corrugation  would  form  in  the  long,  flat 
section  of  the  meniscus  line  on  the  left-hand  of  figure  46(a)  and  a 
corrugation  would  slip  off  the  right-hand  side  leaving  a  pattern  of  the  kind 
shown  in  figure  46(b).  Figure  46(c)  shows  an  intermediate  stage  where  a  new 
corrugation  has  formed  on  the  right,  and  the  corrugation  at  the  extreme  left 
is  in  the  process  of  'slipping  away'.  Unfortunately  I  observed  this  flow  only 


An  unsteady  (V»)  flow  «t  B  *  *1.7,  s  •  2.25,  h  -  8.12  mn,  v  »  0.730  St.  (a)  Five  corrugations  accumulated 
near  the  right-hand  wallj  I  hi  flva  corrugations  accumulated  naar  the  left-hand  walls  Ic)  a  phase  when  at* 
corrugations  war*  present. 
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for  a  20  minute  period  and,  on  reflection,  it  is  feasible  that  eventually  it 
might  have  evolved  into  another  kind  of  motion. 

A  further  small  reduction  in  flow  rate  brought  about  another  structural 
change  with  the  flow  settling  into  a  steady  (VF(5))  motion.  The  evolution  to 
this  state  is  depicted  in  figure  47:  first  a  corrugation  'slipped  away'  at 
the  right,  as  shown  in  figure  47(a)  and  (b);  at  this  stage  it  seemed  (cf. 
figure  47(c))  as  though  the  vertical-face  structure  would  collapse  in  favour 
of  a  lower-face  flow;  but,  as  seen  in  figure  47(d),  a  fifth  corrugation 
developed  and  eventually  the  pattern  equilibrated  itself  forming  the  steady 
motion  seen  in  (e).  (A  side  view  of  this  particular  flow  has  already  been 
shown  in  figure  4(d).)  An  interesting  feature  of  the  steady  (VF(5))  patterns 
was  that  the  meniscus  line  at  the  troughs  of  the  corrugations  passed  round  the 
lower  corner  onto  the  underside  of  the  plate,  giving  the  appearance  of  a 
buttressing  of  the  rear  face  of  the  liquid  sheet.  The  branch  of  steady 
( VF ( 5 ) )  flows  was  stable  to  further  reductions  of  the  flow  rate  until  the 
conditions  R  =  21.5,  S  =  2.33,  h  =  7,98  mm,  v  0.754  St.  were  reached.  The 
steady  (VF(5))  flow  was  observed  for  about  15  minutes  at  this  setting,  and 
appeared  to  be  stable,  but  a  further  small  reduction  in  flow  rate  brought 
about  a  loss  of  stability  of  the  flow  and  the  eventual  formation  of  a  steady 
(LF2)  flow.  The  evolution  to  the  (LF2)  flow  was  a  slow  process,  taking  over  1 
hour  to  complete.  Shortly  after  the  flow  rate  was  changed,  a  corrugation 
'slipped  off'  the  left-hand  side  of  the  pattern,  leaving  a  flow  of  the  form 
shown  in  figure  48(a),  with  five  'lobes'  of  air  protruding  onto  the  end 
face.  Then  each  of  the  end  lobes  disappeared  and  the  three  remaining  air 
pockets  spread  sideways  (see  figure  48(b)).  Eventually  the  extreme  right-hand 
lobe  disappeared,  leaving  the  pattern  shown  in  figure  48(c)  and  then  the  flow 
quickly  reverted  to  the  lower-face  flow  shown  in  figure  48(d). 
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riOJPt  40.  h  ssqoanca  shoving  tha  transition  froa  a  (Vf)  flow  to  an  (Lr2)  flow  at  0  ■  20.8,  8  *  2.37,  h  -  7.92  wm,  V  ■ 
0.758  ft .  (a)  fits  lobaa  prssant t  (b)  photograph  15  alnataa  aftar  (a)  showing  thraa  lobasr  ( o )  photograph 
47  Blnutaa  aftar  (b)  showing  two  lobosi  (d)  photograph  1  alnuts  aftar  (c)  showing  tho  (LT2)  flow. 
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5.  COMMENTS 


(i)  No  attempt  has  been  made  here  to  provide  a  mathematical  discription 
of  the  observed  phenomena  and,  in  general,  this  would  appear  to  be  a  daunting 
task.  It  may,  however,  be  possible  to  give  a  mathematical  account  of  some  of 
the  phenomena  by  studying  a  related  abstraction  of  the  more  complete  model 
outlined  in  §3.  For  example,  the  main  characteristics  of  the  steady  (VF(0)) 
flow,  in  the  zone  near  the  centre  of  the  channel,  should  be  well  represented 
by  a  two-dimensional  abstraction  of  the  more  complete  model.  In  particular, 
the  apparent  uniqueness  of  the  (VF(0))  flow  at  the  larger  flow  rates  of  the 
experiment  suggests  that  such  a  model  could  have  considerable  relevance  to  the 
observed  flow.  Moreover,  should  such  a  flow  exist,  it  would  be  interesting  to 
examine  its  stability  properties  with  decreasing  flow  rate,  to  look  for 
bifurcation  to  a  spatially-periodic  flow  similar  to  the  ’corrugated'  flows 
described  above. 

(ii)  The  experiments  suggested  that,  at  sufficiently  small  flow  rates, 
all  the  observed  motions  were  time  dependent.  It  was  only  at  larger  flow 
rates  that  steady  motions  were  observed.  This  suggests  the  theoretical 
question  of  whether  there  are  any  steady  flows  in  the  limit  that  obtains  as 
the  flow  rate  is  reduced  to  zero. 

(iii)  Many  unsteady  flows  were  observed,  of  which  a  large  proportion  were 
nearly  periodic  in  time. 

Some  of  these  were  the  'drip  flows'  found  at  the  smaller  flow  rates. 

These  flows  had  the  appearance  of  a  basic  periodic  motion,  modulated  over  a 
larger  time  scale  because  the  period  for  droplet  formation  differed  slightly 
from  site  to  site.  For  the  'periodic'  part  of  this  motion  it  seemed 
empirically  as  though  the  relative  phase  of  formation  of  the  droplets  could  be 
fixed  in  any  desired  order  and,  moreover,  it  seemed  as  though  the  phase  at  one 
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drip  site  could  take  any  value  within  some  interval,  relative  to  that  at 
another  site. 

Note  that  a  mathematical  representation  of  the  dripping  process  involves 
the  consideration  not  only  of  domains  which  are  time  dependent,  but  also  of 
domains  whose  connectivity  is  a  function  of  time. 

Many  of  the  'corrugated'  (VF)  flows  had  a  nearly-periodic  structure  in 
time.  The  time-dependent  properties  of  these  flows  seemed  to  be  intimately 
linked  to  the  presence  of  the  side  walls  of  the  reservoir,  since  any 
asymmetries  in  the  pattern  of  corrugations  meant  an  unbalanced  distribution  of 
surface-tension  forces,  resulting  in  one  of  the  corrugations  moving  to  a  side 
wall  and  losing  its  identity.  The  only  corrugated  flow  that  appeared  to  be 
steady  was  the  one  in  which  the  meniscus  line  extended  to  the  underside  of  the 
plate  at  the  trough  of  each  of  the  corrugations. 

(iv)  No  attempt  was  made  in  this  work  to  study  the  dependence  of  the 
flows  on  changes  of  the  geometric  parameters  a,  b,  d  and  a.  It  was, 
however,  evident  in  the  course  of  the  experiments  that  the  flow  structures 
were  strongly  dependent  on  d,  the  distance  of  the  plate  above  the  free- 
surface  level  in  the  reservoir. 

(v)  At  sufficiently  large  flux  rates  there  is  a  potential  long-wave 
instability  of  the  free  surface  in  the  main  channel.  Theoretically  (e.g.  see 
Yih  1969,  p.  502)  this  is  possible  at  flow  rates  such  that  R  >  -g-  cot  a,  i.e. 
when  R  exceeded  approximately  11.3  for  the  channel  slope  used  in  the  present 
experiments.  Thus,  any  unsteadiness  in  the  approaching  flow  from  this  source 
would  have  influenced  only  the  (LF2)  flows  and  the  vertical-face  flows. 
Moreover,  from  monitoring  the  depth  of  the  stream,  the  influence  of  this 
instability  was  not  noticeably  evident  above  the  level  of  the  depth  variations 
observed  at  much  smaller  flow  rates. 


(vi)  Several  of  the  observed  unsteady  flows  had  a  temporal  structure  much 


less  regular  than  the  nearly-periodic  flows  mentioned  above  in  (iii).  They 
all  had  a  certain  recurrence  property  associated  with  them,  but  the  period  of 
the  'cycle*  could  be  dramatically  different  from  cycle  to  cycle,  and  the 
motion  realized  during  one  cycle  often  differed  significantly  from  that  which 
evolved  during  the  next  cycle.  Nevertheless,  the  flow  would  eventually  return 
to  a  structure  that  appeared  to  be  only  marginally  unstable.  In  some  cases, 
(eg.  the  (T),  (C+T),  ...  , flows)  the  recurring  structure  was  one  for  which,  at 
slightly  different  parameter  settings,  there  was  a  steady,  stable  flow.  Thus, 
the  question  arises  as  to  whether  or  not  the  related  mathematical  problem, 
subject  to  steady  flux  conditions,  has  unsteady  solutions  of  the  kind  observed 
experimentally.  For  it  is  quite  conceivable  that,  near  the  margin  of 
stability  of  a  flow,  temporal  imperfections  could  move  the  flow  parameters 
back  and  forth  between  the  basin  of  attraction  of  one  steady  solution  and  that 
of  another,  resulting  in  a  forced  unsteady  motion,  driven  by  imperfections  in 
the  boundary  data.  Whatever  the  source  was  for  these  unsteady  flows  they 
would  appear  to  be  of  considerable  importance,  for  such  unsteadiness  could 
well  arise  in  other  flow  situations  where  there  is  a  large  degree  of 
multiplicity.  Note  that  several  of  the  observed  unsteady  flows  (e.g.  the 
(TCT) ,  ( TCTC ) ,  ...  , flows)  apparently  had  no  associated  steady  'base' 
patterns;  and  one  of  these  flows,  namely  the  (TCT)  flow,  lost  stability  to  yet 
another  unsteady  motion. 
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